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CHAPTER 1: INTRODUCriON 
1.1 Preliminary Remarks 

Metal fatigue under stress and thermal cycling is expected 
to be a principal mode of failure in such engine components as turbine 
blades and disks, and combustor liners. But fatigue design factors are 
subject to considerable uncertainty. For example, enormous scatter is 
observed in cycles to failure fatigue test data with coefficients of 
variation (standard deviation/mean) ranging typically from 20 to 70%. 
Furthermore, scatter exists in operating enviornment data, and uncertain 
ties are present in the models used to predict stresses. Therefore life 
predictions, which rely on the fatigue models and the data for such 
models, are also subject to unde'^tainty . It is suggested that the 
appropriate mathematical model to describe fatigue design factors is a 
probabilistic one rather than a deterministic one. Uncertainty in en- 
vironmentc and in fatigue resistance imply uncertainty in fatigue life 
predictions. This uncertainty can be analyzed rationally only using 
probability theory. 

A reliability approach to high temperature fatigue using 
probabilistic design theory has at least the promise if not the guaran- 
tee of producing better engineered design, i.e., components which are 
more safe, reliable, and cost effective, relative to a deterministic 
approach. Typically conventional design procedures tend to be 
conservative and produce inconsistent levels of risk in components of a 


system. The payoff for an improved design criteria, e.g., for engine 
components, would be a savings in weight. 

1 .2 General Comments on Mechanical Reliability 

The various uncertainties which exist in a mechanical de- 
sign problem can be divided into statistical and nonstatistical (or 
professional) uncertainties as follows 

1. Statistical Uncertainty (Data generally available or easy to 
obtain) 

a. In basic material behavior, e.g., the scatter observed 
in basic S-N fatigue data. 

b. In the estimates of the design parameters from data, 
e.g., the estimate of the fatigue strength coefficient is a 
random variable having a significantly large variance when only 
a small sample of fatigue data is available. 

c. In the mechanical environment (load, high temperature, 
corrosion, etc. which affect fatigue life) 

2. Nonstatistical Uncertainty 

(associated w’th assumptions made in the analysis) 

a. In the theoretical model used to describe fatigue 
strength, e.g., linear damage accumulation rules. 

b. Introduced by the procedure used to compute nominal 
stresses and temperatures in a component, i.e., assumptions 
made in the computer analysis routines. 

c. Introduced by the nradels used to calculate fatigue 
strains at the critical points. 
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d. In strength due to size effects, processing and machin- 
ing operations, assembly operations, etc. 

The above illustrations suggest that we are dealing with a 
highly unpredictable physical phenomena. Some questions immediately 
come to mind. 

1. Do present analysis procedures provide results which are too 
conservative? Are components overdesigned? 

2. Are analyses telling us to provide too little material, or to 
recomnend improper details to insure adequate protection for fatigue? 

3. In general, are the procedures which are currently being used 
really producing high quality designs? 

4. Can we quantify changes in risk (or probability of failure, or 
percent items which fail) due to modifications in the design? 

There is a need to take a hard look at the high temperature, 
low cycle fatigue design process from a probabilistic viewpoint to attempt 
to answer these questions and thereby formulate design strateaies which 
ould reduce engine weight and provide real cost benefits. 

Because of significant uncertainties in fatigue design fac- 
tors, a probabilistic-statistical approach seems particularly relevant. 

All design factors are treated as random variables. Recent develop- 
ments in probabilistic design theory can be utilized to predict distri- 
butions of fatigue life, or to establish design rules. 

Commonly slated reasons for using a probabilistic approach 
to asign include the following: 
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1. It is argued that reliability (or its complement, risk or 
probability of failure) is the most meaningful index of structural 
performance. 

2. The effect on risk of making a design modification can be 
quantified. 

3. A mechanism is provided for explicitly accounting for avail- 
able statistical data on design factors, e.g., in cycles to failure 
fatigue data. 

4. Factors which have nonstatistical uncertainty (due to assump- 
tions made in analysis) can be treated as random variables. Their ef- 
fect on design can be quantified and their relative importance assessed. 
This information can provide guidance for decision makers regarding 
which elements of the problem require further scrutiny. A rational 
basis for decisions regarding research programs is available. 

5. All components can be designed to a balanced level of risk, 
thereby producing a more efficient system. 

6. Probability based information on mechanical performance can be 
used to develop rational policies towards pricing, warranties, spare 
parts requirements , etc. 

In summary, probabilistic design theory provides an improved 
engineering representation of reality. 
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1 .3 An Illustration of the Conservatism in Conventional Design 
Procedures 

Design procedures and codes traditionally rely on a factor 
of safety (applied to material strength) to insure acceptable reliability. 
Such procedures have oerformed well, but it is generally thought i.hat 
regui rements are typically overly conservative and that they produce 
designs having an inconsistent level of risk among the components. At- 
tempts to estimate materials savings which could be realized by prob- 
ability based design criteria have produced figures of about 10%. Thus, 
application of reliability methods to venicle power plant iesign promises 
a payoff in weight savings at the same reliability levels. 

As an example, consider the simple design problem illus- 
trated in Fig. 1.1. The bar is subjected to a harmonic axial load Q(t). 
Fatigue is considered the principal mode of failure. The service life 
is given as = 500,000 cycles. 

Because of analysis procedures, the amplitude Q is con- 
sidered to be random. Statistics on Q are available. The mean value 
is uq = 16 kips (71.2 kN) and the coefficient of variation Cg = 0.10 . 
The distribution of Q is lognormal. As a design value, a point on the 
safe (upper) side of the distribution is chosen as the median plus 
three standard deviations on a log basis. The calculations are sum- 
marized in Table 1.1- See also Appendix I. 

Data is available on the fatigue behavior of the material as 
shown in Fig. 1.^- A basic linear model is assumed (See Chapt. 3 and 4) 
on a log basis, and the least squares curve established. The model for 
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Fig. 1.1 


A DEMONSTRATION OF HOW CONVENTIONAL DESIGN 
CAN PRODUCE CONSERVATIVE RESULTS 



Q(t) 



Qq * 21.5 Kips 
/ (95.7 kN) 

h/ Uq * 16 KIPS 
(71.2 kN) 
Cq = 0.10 

TIME 



CYCLES TO failure. N 


cycles to failure is N = where the median and coefficient of varia- 
tion of K, denoted as K and Cj^ are given in the figure. 

But the designer chooses a curve on the safe (lower) side 
of the distribution as the least squares line minus three standard 
deviations on a log basis. This curve is defined by NS”’ = . 

The conventional approach commonly will require a design 
life higher than N^. A factor of 2 is used here. The design require- 
ment using conventional methods is given in Table i.l. 

Because statistics are available, the risk inherent in this 

design can be estimated as shown in Table 1.1. Definitions and details 

are provided in Chapter 7 and Appendix I. The estimated risk of 
£ 

1.92 X 10 is generally considered to be overly conservative for mech- 
anical and structural components where the consequences of failure are 

-3 

not disasterous. A target value of the probability of failure » 10 
is accepted as being more reasonable. 

_3 

Using a value of = 10 as the basic requirement, a 
probability based design produces a smaller component. The estimated 
weight savings of 15% by using a reliability approach compares with an 
often quoted estimate of 10% for the level of excess material required 
by present design codes. 

In summary, it was intended that this example illustrate 
the fact that excessive levels of reliability can be produced by a 
"pile-up" of safety factors applied to each component of the design al- 
gorithm. It was not the intention to imply that all deterministic 
design codes and procedures are overly conservative. 
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TABLE 1.1 


Calculations: Conventional vs. Probabilistic OesiQi^ 


CONVENTIONAL (DETERMINISTIC) DESIGN 

•Establish Design Values (see also Appendix I for mathematics of 

lognormal vanates) 


LnQ = inQ + 3o, 


hQ 






= 21.5 kip 
= 95.7 kN 



= 2nK - 
"znK = 


K = 1 .25 X 10^® (ksi units) C 
0 

= 5.20 X 10^^ (MPa units) 


• Design Life 

N = 2N = 100,000 cycles 
0 s 

• Design Stress 

$0 = Qq/A = 21.5/A ksi 

• Design Equation 

Cycles to failure, N = K^S^ ^ '^0 

• Solution 2 2 . 

Area, A ^ 0.78 in (5.03 cm ) 

Diameter D ^ 1 .00 in (2.54 cm) 

RISK IMPLIED BY THIS DESIGN 
• Safety Index 
3 = ^n(N/N3) 

'^znN 

N = K/(Q/ ■ )'^ 

-„N ' 

^ = 4.62 


•Probability of Failure 

(-.'!/ 


p^ = 1.92 X 10 


= 0.40 
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. . . TABLE ^ . 1 Continued 


REDESIGN USING RELIABILITY METHOD 

•Target = 10 (reasonable level) 

Target o = 3.09 

•Solution (using above equations) 

Area, A = 0.66 in’ 

Diameter D = 0.92 in 

PERCENT WEIGHT SAVINGS BY USING RELIABILITY METHOD 
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1 .4 Goals of This Stud y 

The ultimate goal of this study is to recommend methods of 
reliability analysis for the development of probability based design 
criteria for fatigue design in general and high temperature fatigue in 
particular. It is intended that this study provide practical information 
for engine designers. 

Specific goals include the following. 

1. To develop methods of providing statistical summaries of 
data for design purposes, i.e., 

a. Design values and/or design curves 

b. Appropriate statistics which are required for a 
reliability analysis. 

2. To recommend methods of probabilistic design for a 
complete reliability analysis of a component subjected to 
high temperature, i.e., the problem where several random 
design factors appear in an equation and where it is re- 
quired to evaluate the risk, 

3. To present the results of this study in a format for easy 
implementation for engine designers . 
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1.5 What is Contained in This Report 


This report summarizes the results of the first year of a compre- 
hensive study on fatigue/creep reliability supported by the NASA/Lewis 
Research Center, It is intended that this report provide practical 
information useful for designers. 

Chapter 2 summarizes methods available for statistical analysis 
of ata on a single variable. Various schemes for establishing a design 
value are presented. 

Chapters 3, 4 and 5 summarize procedures for providing statistical 
summaries of S-N fatigue data. Linear model analysis is presented and 
applied to both homo- and heteroscedastic data. Preliminary consider- 
ations of data analysis using the general strain-life model are included. 

Attempts are made in Chapter 6 to present a coordinated overview of 
the three basic approaches to formulating a fatigue design equation. 
Performance of Miner's rule and strainrange partitioning are described 
by statistical summaries. 

Chapter 7 probides a summary of available reliability methods for 
fatigue design. Two examples are presented in Chapter 8. One is a 
strainrange partitioning example. The other is of a local strain analysis 
model . 

This report is considered preliminary. Work is continuing in all 
areas, and future reports will provide more complete information. 
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1 . 6 List of Genera l Refere n ces for Reliability and Probabilistic Design 
The following references include some of the more important 
sources of information in elementary applied proDability theory and 
statistics, probabilistic design and fatigue reliability. The list 
is by no means exhaustive and does not include many excellent texts 
in specialized areas of reliability engineering, quality control, 

random process theory, etc. 

for basic probability and statistics 

Ang, A. H.-S, and Tang, W. , Probability Concepts in Engineering 

Planning and Design , Wiley, 1^/5. Excellent text tor engineer- 
ing problems and statistics, but oriented to civil engineering. 

Benjamin, J. R., and Cornell, C. A., Probability, Statistics and 
Decision for Civil Engineers , McGraw-Hi 1 1 , An excellent 

reference for engineers oi^ all types, uses modern techniques 
with many examoles, probably one of the best general references 
available for engineers to date. 

Bowker, A. H., and Lieberman, G. J., Engineering Statistics , 

Prentice-Hall, 1972. An excellent reference, strong on test- 
ing hypotheses and quality control. 

Hines, W. H., and Montgomery, D. C., Probability and Statistics in 
Engineering and Management Science , 2nd Ed., Wiley, 1980 

Lipson, C., and Sheth, M. J., Statistical Design and Analysis of 
Engineering Experiments , Mc6ra\w-Hill , 1973. Treats problems 
from mechanical engineering design theory. 

Meyer, P. L., Introductory Probability and Statistical Applications , 
Addi son-Wesley, 157^1 An excellent introduction to modern 
probability theory and mathematical statistics. 

Mood. A. , and Gr^ybill, F. A., Introduction to the Theory of 
S_tdtisti«, McGrew-Hill, 1963. Intermediate statistics, and 
a~TTttl e difficult for a beginner. 

Natrella, M. G., Experimental Statistics , Oeot. of Commerce, 1963. 

An excellent statistics "cookbook" full of examples on how-to- 
do-it in testing hypotheses and many other tests. 
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for probabilistic design 


Ang, A. H.-S., and Cornell, C.A., "Reliability Basis for Structural 
Safety," Journal of the Structural Division , ASCE, Sept. 1974. 

Cornell, Ca.A., "A Probability Based Structural Code," ACI Journal, 

Oec. 1969. 

Dialog . Second International Workshop on Code Formats, Mexico City, 

Jan. 1976, Danmarks Ingeniorakademi , Building 373, 2800 Lyngby, 
Denmark. 

Ellingwood, B., et al., "Development of a Probability Based Load 
Criterion for American National Standard ASS," NBS Special 
Publication 577, June 1980. 

First Order Reliability Concepts for Design Codes , Bulletin D' Infor- 
mation 112, Comite Europeen du Beton, Munich, July 1976. 

Freudenthal, A.M., Garrelts, J.M., and Shinozuka, M., "The Analysis 
of Structural Safety", Journal of the Structural Division, ASCE. 
Vol.92, No STl, Feb. 1966: 

Galambos, T.V., Ravindra, M.K., (and others),... a series of eight 

papers on Load and Resistance Factor Design (LRFD) in the Journal 
of the Structural Division . Vol. 104, No. ST9, Sept. 1978. 

Hasofer, A.M. , and Lind, N.C., "An Exact and Invariant First Order 

Reliability Format," Journal of the Engineering Mechanics Division, 
ASCE, February 1974. 

Haugen, E.B., Probabilistic Mechanical Design . Wiley, 1980. 

Haugen, E.B., and Wirsching, P.H., "Probabilistic Design," a five part 
series in Machine Design , starting April 17, 1975. 

Kececioglu, D.B., and Cormier, D., "Designing a Specified Reliability 
into a Component," Proceedings of the Third Reliability and Main- 
tainability Conference, Washington, D.C , 1964. 

Probabilistic Mechanics and Structural Reliability , (ed. A. H.-S. Ang 
and M. Shinozuka), ASCE Tucson, January 1979. 

Probabilistic Methods in Structural Engineering (ed. M. Shinozuka and 
J.T.P. Yao), ASCE, 1981. 

Rationalisation of Safety and Ser viceability Factors in Structural 

Codes, Report 63, CIRIA, Construction Industry Research and Infor- 
mation Association, 6 Storey's Gate, London SWIP 3AU. 
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for fatigue reliability (many other references cited in text) 

Ang, A. H.-S., and Munse, iJ.H., "Practical Reliability Basis for 

Structural Fatigue," ASCE National Structural Engineering Confer- 
ence, Preprint 2494, April 1975. 

Ellingwood, B., "Probabilistic Assessment of Low Cycle Fatigue Behavior 
of Structural Welds," Journal of Pressure Vessel Technology , ASME, 
Paper 75-PVP-29, February 1976. 

Fatigue Reliability; A State of the Art Review, a four part series 

1. "Fatigue Reliability: Introduction" 

2. "Fatigue Reliability: Quality Assurance and Maintainability" 

3. "Fatigue Relibility: Variable Amplitude Loading" 

4. "Fatigue Reliability: Development of Criteria for Design" 

Journal of the Structural Division, ASCE, Vol . 108, No. STl, Jan. 
1982. 

Fong, J.T., "Uncertainties in Fatigue Life Prediction and a Rational 
Definition of Safety Factors", Nuclear Engineering and Design , 

Vol. 51 , 1978. 

Freudenthal, A.M., and Gumbel, E.J., "On the Statistical Interpretation 
of Fatigue Tests," Proceedings of the Royal Society of London, 

Series A, Vol. 216, 1953, pp. 309-322. 

Freudenthal, A.M. and Gumbel, E.J., "Physical and Statistical Aspects 
of Fatigue," Advances in Applied Mechanics , Vol. 4, 1956. 

Kececioglu, D.B., and Chester, L.B., "Combined Axial Stress Fatigue 

Reliability for AISI 4130 and 4340 Steels," ASME Paper 75-WA/DE-117, 
1975. 

Little, R.E., and Jebe, E.H., Statistical Design of Fatigue Experiments , 
Applied Science, 1975. 

Whittaker, I.C., Besuner, "A Reliability Analysis Approach to Fatigue 
Life Variability of Aircraft Structures," Wright-Patterson Air Force 
Base Technical Report AFML-TR-69-65, A;ril 1969. 

Wirsching, P.H., "Fatigue Reliability of Welded Joints in Offshore 
Structures," Proceedings of the Offshore Technology Conference OTC 
3380, 1979, also published in the International Journal of Fatigue , 
April 1980. 

Yang, J.N., "Statistical Approach to Fatigue and Fracture Including 
Maintenance Procedures" Fracture Mechanics , U. of Virginia Press, 
1980. 

Yao, J.T.P., "Fatigue Reliability and Design", Journal of the Structural 
Division , ASCE, Vol. 100, No. ST9, September 1974. 
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Chapter 2 METHODS OF DATA ANALYSIS OF A SINGLE VARIABLE 


2.1 Preliminary Remarks 

Given a random set of observations on a design factor, it is necessary uo 
interpret the data in order to make a design decision. For example the random 
sample of Table 2.1 is cycles to failure fatigue data for tests at a single 
stress level. In order to insure an acceptably low level of risk (of proba- 
bility of failure) it may be necessary to specify either a "design value" on 
the safe (lower) side of the distribution or the statistics of the variable, 

depending upon the design strategy used. 

This chapter summarizes various statistical tools which can be employed 
to provide designers with information that they need to make decisions on a 
single variable. Later chapters deal with analysis of SN data. The focus of 
this report is on fatigue, but techniques described herein have wide appreciation. 

2.2 Mathematical Tools for Probability Estimates 

Because observed cycles to failure has significant scatter, it is suggested 
that cycles to failure, denoted as N, be treated as a random variable. There- 
fore a probability density function (pdf) fjyj(n) is defined such that 

P(n^ < N < ng) = jfj^(n)dn (2-1) 

"a 

where P( • ) denotes the probability of the event in parentheses. 

The cumulative distribution function (cdf) Fj^(n) is defined as 

n 

F^(n) = P(N < n) = jf,^(x)dx (2.2) 

— CO 

where n denotes a specific value of the random variable N. 
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Table 2.1 


Cycles to Failure of Specimens Tested at the Same Stress Level [Ref: Evans (13)] 


3 

Cycles to Failure (10 cycles) 
15.4 
22.2 

17.3 

23.6 

14.4 

12.3 

16.5 

25.7 

17.5 

27.0 

20.5 

21.3 

14.0 

27.0 

23.0 

23.6 

16.0 

14.7 

10.3 
16.0 

13.0 

31 .0 


Sample Size, n 
Sample Mean, N 

Sample Standard Deviation, Sj^^ 
Sample Coefficient of Variation 


22 

19.2 (10^ cycles) 
5.59 (10^ cycles) 
29. ir. 
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An example of a pdf and corresponding cdf is given in Figure 2.1. Note 
that all of the information about probabilities of N is contained in the pdf 

and cdf. 

Tn mean and standard deviation are defined as. 


Pm = I nf^(n)dn 


-T58- 


"n = 


J 


(n - u^) ff^(n)dn 


(2.3) 

(2.4) 


The mean is a measure of the central tendency and the standard deviation is 
an index of the degree of variability. Another measure of central tendency is 
the median, N, defined as the 50« point of the distribution 


Fj^(N) = P(N < N) = 0.50 


(2.5) 


The coefficient of variation is defined as 


This term, commonly used in probabilistic design literature, is a non-dimensional 
measure of variability. 

Ultimately we wish to make statments regarding risk (or probability of 
failure) which we can do if we have F^^(n) and/or f|^(n) . But at this stage all 
we have is the random sample N., i = l,k of Table 2.1. The remainder of this 
chapter deals with the problem of statistics, i.e. how we take the random sample 
and make probability statements for design purposes. The example used herein 
is of cycles to failure data, but the analysis presented applies to any random 

variable. 
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2.3 Basic Statistics 


The sample mean N is a measure of the central tendency 


N = 


1 - 

i=1 


( 2 . 6 ) 


where k = sample size. H is an estimate of 

The sample standard deviation is a measure of the dispersion or scatter 


in the data 


Sm = 




K 


N)‘ 


1/2 


(2.7) 


is an estimate of 

For the data of Table 2.1, N = 19.2 and * 5.59. The sample coefficient 
of variation is = Sj^/N = 0.291. 

This value of C|^ is typical of fatigue data at relatively low lives as 
illustrated in Table 2.2. The values listed give evidence of the relatively 
large scatter which exists in fatigue strength data. For example, yield and 
tensile strengths for a wide variety of materials are typically less than 10% 
and usually about 5%. 

The data of Table 2.2 also suggest that scatter is greater at lower stress 
levels and longer lives. This seems to be a general rule for smooth specimen 
data. 

The emlirical cumulative distribution function, an estimate of Fj^j(n), can 
be established as demonstrated in Table 2.3 (the data from Table 2.1). A 
smooth curve through these points is an estimate of Fj^(n). Examples of how 
the empirical cdf is used are given in Section 2.5. 
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Typical Values of the Coefficient of Variation in Cycles to Failure Data 
Observed by Various Investigators at Different Stress Levels 
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Numbers in parentheses refer to references at the end of each chapter. 


Table 2.3 


The Emprical Cumulative Distribution Function for the Data of Table 2.1 

(n = 22) 


Order Number, i j 

Ordered Data, 

* i - 1/ 2 

Empirical cdf > ^ 

1 

10.3 

.023 

2 1 

i 

12.3 : 

.068 

3 1 

13.0 

.114 

4 

14.0 

.159 

5 

14.4 

.205 

6 

14.7 

.250 

7 

15.4 

.295 

8 

16.0 

.341 

9 

16.0 

.386 

10 

16.5 

.432 

11 

17.3 1 

.477 

12 

17.5 

.523 

13 

20.5 

.568 

14 

21.3 

.614 

15 

i 22.2 

.659 

16 

1 23.0 

.705 

17 

23.6 

.750 

18 

1 23.6 

.795 

19 

1 

25.7 

.841 

20 

27.0 

.886 

21 

27.0 

.932 

22 

31.0 

.977 


I i 


*Many other forms which give similar results have been suggested, but this 
one (sometimes called the "Hazen formula") seems to generally perform well 
and is widely used by engineers. 
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- • Stati 5ti c al _ M odels ijs ed i n Fatigue 

In order to make design decisions based on a set of observations of a design 

factor, it is necessary to describe the distribution of that factor. In that 
regard, statistical models are usually employed. The random variable N denoting 

cycles to failure is often described with a two parameter Ueiruli or 
lognormal model. A summary of these models as well as the normal and the 
three parameter Weibull is provided in Tables 2.4 through 2.7. The normal is 
included in this discussion primarily for reference. The three parameter Weibull 
has been used as a model for fatigue data; practical and analytical difficulties 
associated with the use of this distribution are summarized later. 

The lognormal and two parameter Weibull are most commonly used to describe 
N. Use of the lognormal distribution has been based primarily on arguments of 
mathematical expediency. However physical arguments favor the Weibull for 
most material strength variables because it is an asymptotic distribution of 
minima of a sample (/). If failure of a structural element is precipitated by 
failure of the first of a large number of sub-elements, then the Weibull is 
1 ikely a "good" model . 

Moreover, it has been pointed out by Gumbel (8) and demonstrated in Figure 
2.2a that the hazard function decreases for large M which violates our physical 
understanding of progressive deterioration resulting from the fatigue process. 

Also mote from Figure 2.2b that for the same statistics, the Weibull gives 
larger probabilities in the left tail. Use of the Weibull should produce 
conservative designs relative to the lognormal. 

Nevertheless, the lognormal is often much easier to use. Methods of 
linear model analysis commonly used on SN data rely on a lognormal assumption 
for N. P“*obabil istic design procedures in a lognormal format have been 
developed. Furthermore, this author has found, more often than not, the log- 
normal provides a better fit than the Weibull to real fatigue data. 
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TABLE 2.4 

A SUMMARY OF THE NOR'^AL DISTRIBUTION 


Probability Density Function 


TM(n) = exp 


r_ 

I 

L 




J 


Statistical parameters 

U|,j = mean value of N 

Onj = standard deviation of N 


Distribution Function; Probability Calculations 


n^ - u 


= P(N < n^) ° ^ 




L J 


where f ( • ) = standard normal distribution function 
(definition, e.g. Reference 14) 
n^ = any specific value of N 


How the Parameters are Estimated from Data 

Data : N = (N.| , . . . N|^) 

compute the sample mean M as an estimate of 

comoute the sample standard deviation Sj^ as the estimate of 
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Table 2.3 

A. SUVMAKY OF THh LOGmQWal 0ISTF1IBUTI0:J* 


‘1 1 -; Said to nave a lognormal distribution if X = In N (or '( = loo Ti) has a 

■ lO 

normal distribution. 


-Probability Oensi ty Function (base e) 


f.(n) 


N’"' 


!" (In n - -M )“ 

exn ! - ^ 


cCT 


T 


Statistical parameters 


= mean value of X 
a = standard deviation of X 

X 


•Distribution ^unction Probability Calculations 


F.|(n^) = P(» < n^) = t| "o ■ “x'' 


or 






V. 


J 


'vhere ?(*) = standard normal distribution function 


*See Appendix 1 for a comprehensive summary of the lognormal distribution 
and its application in probabilistic design. 
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Table 2.5 (continued) 


.Relationships between Parameters and "oments (i.e., mean and standard deviation) 


Base e 

-X " ’"N* + C‘) 

= In N (See note below) 

= In (1 + C‘) . 


= expCuj( + 


- 1) 

C„ =yexp(o^2) _ , 


Base 10 




aj = 0.434 log^^d + C‘) 

. lol-x * l/2(aj/.43<!)t 

/ (aJ/.A34) 

Cm = ;io ^ - 1 


How the Parameters are Estimated from Data 

Data N = (N^.n^.-.Ni^) 

Let X^. = In N^. (or loq.jQN^. ) 

Let sample mean X be tne estimate of 

Let sample standard deviation s be tne estimate of a 

X X 

Note: A useful relationship between the median N and the mean U|y 

can be derived as ^ 
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Table 2.6 

A SUMMARY OF THE Tl>iO PARAMETER WEIBULL DISTRIBUTION 


Distribution Function; Fj^j(n) - P(N < n) 

Fj>^(n) = I - exp l_- (j) j 

Statistical Parameters 

1 = shape parameter 

8 = scale parameter 


Probability Density Function 

fj^(n) = (a/’3)(n/s exp |-(f) ! 


Relationship between Parameters and Moments 


u 


N 


C 


N 


3r(1/a + 1) 



where r(.) = the gamma function 

- 1 08 

A useful approximation a ^ ' for 0.02 < Cj^ < 2.00 


How to Estimate the Parameters (Method of Moment Estimators) 

Data N = (N^ , N 2 , . • • \) 

Compute sample mean R and sample standard deviation of Sj^^ 

Let a ={s,^/N)'^ ®n estimate of a and 3 = N/r(1,/ot + 1) be an 

estimate of 3 
(Method of Moments) 
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Table 2.7 


A SUVrARY OF THE THkEE PARa^^ETER WEIBIlL 


Distribution Function; F^,(nj = < n) 


I ! 




F|,j(n) = 1 - exp j 


a 


J 


Probability Density Function 


fr,,m) = fl^^j exn I - 




L 


Relationshio Between Parameters and *^oments 


Uj^ = Y + sr(i + l/a) 

= 3'[r(l + 2/a) - r^(1 + 1/a)] 

How to Estimate the Parameters 
Maximum Likelihood Estimators 
Ref; Mann et al (9) 

Least Squares Estimators 
See Sections 2.5 and 2.6 
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FIGURE 2.2 COMPARISON OF THE LOGNORMAL 



2 . 5 Probability Plots; Rectification of the Distribution Function 

A probability plot can ba usad as a basis for making a subjective judgement 
on how well a set of data fit a particular distribution. Following are the 
transformations on the normal, lognormal and VJeibull (2 and 3 parameters) which 
permit their distribution functions to be plotted as straight lines. 

(a) Normal distribution ; the distribution function is. 


F„(n) 



(2.S) 


where is the standard normal cdf, and u and a are the mean and standard 
deviations of the random variable N. 

Inverting, this expression 


a 


(2.9) 


and 


= Isl" - a. 


( 2 . 10 ) 


A linear relationship Y = AX + B has been established '-.'here 


Y = X = n, A = 1/a, B = -u/a 


( 2 . 11 ) 


(b) Lognormal distribution ; the distribution function is 


F^(n) = ^ 


In n - 


( 2 . 12 ) 


where X = In N. Here the statistical parameters are and which are the 
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mean and standard deviations of X = In N. Inverting, 


in n - u, 


■ = ^'^(F) 

(2.13) 

* — ' In N - * ''X ■ 

(2.14) 

.""X, ! -a. i 

A' 



A linear relationship has been developed for and In (n). Here 


Y = f (F), X = In N, A = 1/c^, B = (2.15) 


The mean and coefficient of variation of N is. 


"N “ 


( 2 . 16 ) 


Cjy, = 7exp(cp - 1 

(c) We ibull (2 parameter) distribution ; the distribution function is, 

Fj^(n) = 1 - exp[-(n/S )'■''] (2.17) 

where ot and 3 are the statistical parameters. Inverting, 

-(n/3)^ = In (1 - F). (2.18) 


Taking the log of both sides. 


a(ln n - In 3 ) = In [-In (1 - F)] 


(2.19) 
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Let, 


Y = In [-In (1 - F)] 
X = In N. 


( 2 . 20 ) 


Then, Equation 2.17 becomes linear, 

Y = AX + B (2.21 ) 

where 

A = a, B = -a In 3 (2.22) 


The mean and coefficient of variation of N in terms of a and 3 are 


U„ . sr{- + 1) 


c„ = a 


.925 


(2.23) 


where r(*) is the gamma function. A chart which enables convenient evaluation 
of the gamma function is provided (Figure 2.3). The expression for C,^ is 
approximate, but the error is small for 0.02 < < 2.00. 

(d) Weibull (3 parameter) distribution ; The transformation described above 
for the two parameter Weibull applies directly to the three parameter model 
when the following transformation is made 

Y = ln(N - y) (2.24) 


The parameter v must be known in advance (see discussion below on the least 
squares method). 
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Using the above information for the transformation of each statistical 
model, one can plot F vs. N as a straight line on rectangular paper with the 
appropriate variable representation on both X and Y axes. Incidentally, this 
operation is performed routinely using appropriate probability paper but 
automatic plotting was used herein. 

The empirical distribution function F for the example defined in Tabie 2.2 
was plotted in Figure 2.4 on "lognormal paper". Note that F^. is an estimate 
of the actual cdf, Fj^(n). Thus, because the data tends to plot as a straight 
line in "^igure 2.4, the lognormal may be a "reasonable" model for N. This 

test is of course, entirely subjective. 

Addressed now will be the question of how to determine quantitatively 
which model, e.g. lognormal or Weibull, best fits the data. This issue is 
discussed in Sections 2.6 and 2.7. 

2.6 Comparison Test Using the Least Squares Method 

Using probability plots, comparison of the fit of each model can be made 
analytically by (a) computing the least squares estimators thereby defining 
the least squares line, (b) measuring the amount of deviation of the data from 
the least squares line and (c) choosing that model which has the least deviation 
Using the least square method, a straight line can be fitted through these 
data points (see Sec. 3.2.1). The least squares line is an estimate of 
E(Y|X) = A + uX 

Y = A + BX (2.25) 

where the least squares estimates A and B of A and B are, 

','X Y. - nXY 

B = (2.26) 

Zxi - n(M- 
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Figure 2.4 


Fatigue Data of Table 2.1 Plotted on Lognormal 
Probability Paper 
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A = Y - ex 

and where X and Y are the sample means of X and Y, respecti vely. 

The least squares line is shown for normal, lognormal, and Weibull plots 
of the data of Table 2.2 in Figures 2.5, 2.5 and 2.7. The purpose of such an 
analysis is to provide a test (albeit subjective) for establishing which 
model provides the best fit of the F-N data. The model for which the data 
seems to best plot as a straight line is the one which is the "best fit". 

Such an exercise is not always successful in identifying the best model 
as Figures 2.5, 2.6 and 2.7 illustrate. All three models seem to provide 
a reasonable representation. Therefore, an objective test is necessary. 

A computer program (LESQUE-I) has been developed to analyze F-N data. 
The output includes, (a) determination of the least squares estimators A 
and 8, (b) the parameters of the model from the least squares line and 
(c) the sample correlation coefficient :. The coefficient of correlation 
is used as a measure of goodness fit. The sample correlation coefficient, 
denoted by o , is expressed as. 


n 


i ) 

n 


X. - X 


Y. - 
1 


1 = 1 ^ 




(2.27) 


with -1 < G < 1 and where X and Y are the sample means and 

n _ 


1 


S = - 


n 


(X. 


i = 1 


(2.28) 


Sy has a similar expression. 

c is computed for each model. If all the data fell on a straight line, 

= 1 (or -1). If the data indicates no linear relationship, c = 0. Therefore 

the model having a larger g suggests a better fit. 
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Figure 2.6 

Fatigue Data of Table ‘2.1 Plotted on Lognomal 
Probability Paper 
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The program was used to analyze the data of Table 2.2 for the two 
parameter models* normal, lognormal and Weibull. Upon examination of the results 
in Table 2.8 it is seen that the values of o are relatively close, and the 
lognormal is the winner. A more discriminating index for comparison of these 
data is a correlation parameter defined as 

£ , A . 0^ ' " (2.29) 

A much larger separation exists between these values as shown in the table. 

£ is inversely proportional to o. Thus smaller e indicates a better fit, and 
the lognormal is the best fit of these two parameter models. 

2.7 Least Squares Analysis of the Three Paramete_r 
Ueibun Some Lduorial Comments 

A least squares analysis routine (LESQUE-II) for the three parameter 
Weibull (TPW) was developed independently. Analysis is similar to the two 
parameter case, except that the parameter must be specified before analysis 
can proceed. An iterative scheme is used for obtaining the least squares 
estimator of y, denoted as y. A trial value of y is chosen and r is computed. 
The least squares estimator y is obtained by repeating the process to find 

th6 V3lu6 of V which indk6S 3 iTiaxiiTium. 

An analysis of the data of Table 2.1 was performed. The probability 
plot is given in Figure 2.8. Values of r and £ are given in Table 2.8. 

The three parameter Weibull (TPW) is frequently used as a model for cycles 
to failure (e.g. Fong (16)) and the ASTM has recommended the use of the TPW(17). 
The TPW is attractive because the location parameter y defines a non-zero lower 
bound on the sample space. In theory, such a model seems more realistic than 
the two parameter models which permit values (albeit with small probability) 


down to zero. 
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Table 2.8 


Summary of Results of the Least Squares Comparison Test 
on Data of Table 2.2 


Correlation Coefficient 


Co rr el a t i o n Parame ter 

r: » v"l - 0 ‘ 


1 Parameter Models 



Normal 

j 

.9797 

0.200 

Lognormal 

1 .9876 

0.157 

Wei bull 

.9683 

0.250 

ree Parameter Model 
Wei bull 

.9908 

0.135 




Best fit of the two ;.' 3 rameter models 


Figure 2.8 

Fatigue Data of Table 2.1 Plotted on Wei bull 
Probability Paper (Three Parameter Wei bull) 





There are however some undesirable features of the TPW which may make its 
uSvi impractical in certain cases: 

1. The TPW is a very difficult distribution to use because 

a) A complicated iterative program is required to compute estimates of 
the parameters, e.g. the maximum likelihood estimates (9). 

b) Distributions of the estimators for routine statistical operations 
such as confidence intervals, testing hypotheses, etc. are difficult 
if not impossible to obtain. 

c) Complicated numerical analysis is required (integration or Monte Carlo) 
for any reliability analysis involving any other random design factors. 

2. The results given in Table 2.8 show clearly that the TPW provides a better 
fit to the Table 2.1 data than the three two-parameter models studied, but 
this comes as no surprise. More parameters produce a better fit. Why not 
use a four -parameter model,... or five? 

3. Figure 2.9 show the density function of the fitted TPW relative to the 

data. Values of the location parameter y, and the smallest sample point 
Npj are as shown. Note that y is only slightly below Moreover, 

recall that y is a point below which we are absolutely certain that no 
future values will fall. Examination of these values of y in light of the 
scatter in the small sample suggests that to define y as a lower bound is 
risky to say the least (e.g., the normal model predicts that 3" of the data 
will fall below y) . 

4. Because of the truncation of the distribution at a level which may be too 
high, use of the TPW may result in risk estimates which are too low. Thus 
it is likely that the model may produce unconservative designs. 
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2.P Another Contparisorv Test Based on Fit to the Empirical Distribution Function 
Wirsching and Carlson (10) developsd an objective method for determining 
which of several competing statistical models best describes the data. Their 
test is based on deviations between the empirical and hyoothesized distribution 
functions. 

Let F. denote the empirical cdf, which is determined from the sample and 
is an estimate of the cdf of X. The cdf of the model considered is Gj(x;^) 
where 9 is the parameter as estimated by the data. 

The rationale for the test is that if 6j(x;9) is the sampling distribution 
of X, or a good approximation thereof, and F. will tend to be "close to" 
6.(x;9). The W-statistic, a measure of this closeness, is based upon a form 
similar to the Cramer-von Mises statistic used for goodness-of-fit tests 




1-2 

- I 
n 1.1 


where 


D. . = G.(x.;9) - F. 


(2.30) 


(2.31) 


Deviations associated with the W-statistic are shown in Figure 2.10. 

The value of W. is computed for each of the competing statistical models. 
J 

The model having the smallest W. is considered to be the best fit. The W- 

sJ 

statistic is computed for the normal, locnormal and Weibull using the data of 
Table 2.1. 

W-STATISTIC OF DATA OF TABLE 

0.0591 
0.0491 


Normal 
Lognormal 
Weibul 1 
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0.0525 




This test agrees with the least squares comparison test by ranking the 
lognormal as the best. However the two tests disagree on the normal and Meibull 
It is the opinion of the author that the W-statistic is the more powerful test. 

2.9 The Loqnormal Assumption 

For the most part it will be assumed herein that the lognormal is the 
governing distribution for not only cycles to failure data, but also for the 

other fatigue design factors, see Appendix I for details of method. The 
reasons for the use of the lognormal are; 

1. The lognormal generally has been shown to provide a reasonable description 
for the distribution of a wide variety of design variables. For example, 
upon analysis of cycles to failure data, using methods described above, 
the author has found that the lognormal consistently provides a better 
fit for fatigue data than does the Weibull or normal. 

2. Statistical properties of the lognormal are well defined. All of the 
theory developed for normal, e.g. confidence intervals for the mean, 
apply equally to the lognormal. 

3. The lognormal is easy to use. For example if fatigue life T is a 

multiplicative function of several lognormal random variables, then T 

IS exactly lognormal and it is a simple matter to compute failure proba- 
bi 1 i ties 

4. Reliability formats using the lognormal can easily accomodate design 
variables having relatively large coefficients of variation. Some formats 
rely on small variance assumptions for the design factors. 

5. The lognormal is already widely used in the design profession. For example 
commonly used methods of linear model analysis for characterizing S-N 

fatigue data implicitly assumes that cycles to failure has a loonormal 
distribution (e.c., see Chapte^-s 3 and 4). 
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2.10 The Use of the Tolerance Interval of Establishing a Design Value 
The "tolerance interval" can be used to establish consistent and 


rational design values on the safe side of data. (For a general reference, 
see Natrella (11)). Design values in MIL-HDBK-5, for example, are based 
on tolerance intervals (12). The following discussion describes the 
rationale and demonstrates the operation of establishing the lower toler- 
ance limit as a design value. Unfortunately , while this analysis can 
specify a lower bound to SN data, it falls to provide a mechanism for 
consideration of other factors in the fatigue equation which are subject 
to uncertainty. 

As an example, consider the cycles to failure data of Table 2.1. 

It is assumed that N has a lognormal distribution so that Y=log N has a 
normal distribution. Values of , as well as the sanple mean T and 
sample standard deviation ,s ,of Y are given in Table 2.9. 

Suppose, for example, a decision has been made to establish a design 
value Yq (or Nq= 10 D) as the value below which it is anticipated that no 
more than 1 (a=0.01) of future measurements would be expected to fall. 
Figure 2.11 illustrates relative to the distribution of Y. The value 
of Yp and the corresponding are easily calculated, assuming that Y has 
a normal distribution with mean Y and standard deviation s. 


Design Value for Y : 
In general 



In this example 

Yq = 1.27 - (2.33)(0.127) = 0.9741 


(2.32) 
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Table 2.9 

Cycles to Failure of Specimens Tested at the Same Stress Level 


N. 

1 3 

Cycles to Failure (10 cycles) 


15.4 

1.19 

22.2 

1.35 

17.3 

1.24 

23.6 

1.37 

14.4 

1.16 

12.3 

1.09 

16.5 

1.22 

25.7 

1.41 

17.5 

1.24 

27.0 

1.43 

20.5 

1.31 

21.3 

1.33 

14.0 

1.15 

27.0 

1.43 

23.0 

1.36 

23.6 

1.37 

16 0 

1.20 

14.7 

1.17 

10.3 

1 .01 

16.0 

1.20 

13.0 

1 .11 

31.0 

1 .49 

n=22 

Y=1.27 

s=0.127 
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Figure 2.11 
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Design Value for N : 

In general 

Nq = lO^D (2.33) 

In this example 

Np = 10°-^^^^ = 9.42 (10^ cycles) 

It is estimated that there is only a chance that a specimen 
selected at random would have N cycles to failure less than 9.42 thou- 
cycles. Stated another way, for a large number of specimens we would 

expect about 1% of them to have cycles to failure less than 9.42 thou- 
cycl es . 

But the analysis of Figure 2.11 fails to account for the fact that 
Y and s are themselves random variables. Thus, Y^, as seen from Eq 2.32, 
is a random variable, and Nq from Eq. 2.33 is also a random variable. 

For example, if several labs independently conducted this same test 
on 22 specimens, each would obtain a different value of Nq because o^ 
inherent randomness in the process. Therefore, being a random variable 
has its own distribution function F(Ng) and it makes sense to say, for 
example, that there exists a value of N corresponding to F(Np) = .05. 

This value denoted as N and called a "tolerance limit" is the value above 

a 

which we may predict with 95% confidence that 99% of the population will 
lie. This scheme is coimonly used to establish design values {e.g.. 
Reference 12). 

In general, is the point below which we expect proportion a to 
lie with confidence y- 
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To find a single value Y above which we may predict with confidence y 

Ot 


that a proportion a of the population will lie 


Y = Y - K s 

a a,y 


(2.34) 


where K is found from Table 2.10. Using the above example, n=22, y*0.95, 

a,Y 

a=0.01 and K gg=3.233 from Table 2 . 10 (noting that p=l-a) 

Y = 1.27 - (3.233)(0.127) 

a 

» .8594 (2.35) 


and the corresponding point for N, denoted as for y-.9Sy a*0.01 is 


N » lO^o" » ^qO.8594 

a 


N 


ct * 7.23 thou-cycles 


( 2 . 36 ) 


Consider another example. To establish a "safe life" in turbine discs, 
the following criteria is used for civil engines in the UK 03). The safe life 
is established at minus three standard deviations from the population mean. 
Then there is a 95" confidence that the probability of failure at the safe 
life does not exceed 1 in 750. Note that for the normal distribution, the 
tail area beyond three standard deviations is 0.00135 or approximately 1 in 
750. 

For n=22, 7=0.95 and a=. 00135, the value of K is found from interpo- 
lation of Table 2.10. K gg= 4.20 Using values from the above example, 

= 1 .27 - ( 4.2D',(0.127) 

= 0.7367 (2-37) 


The "safe life" is 

N 

a 


Y 

10 ^^ 


10 


0.7367 


5.45 thou-cycles 


(2.38) 
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TABLE 2J0 FACTORS FOR ONE-SIDED TOLERANCE UMITS FOR NORMAL DISTRIBUTIONS* 


Factors K such that the prob^ility is 7 fhof at leost a proportion P of tho ditfribuHon will bo loss then 
X fCs (or grootor thon X — K$), whoro X and s aro ostimotos of tho moon ond tho stondord 

deviation computed from a somple size of e* 


€ 



y = 0.75 

1 

) 

t 

y - 0.90 

E 

0.75 

0.90 

j 0.95 

0.99 

0.999 

0.75 

0.90 

0-95 

0.99 

0.999 

■ 

1.464 

2.501 

1 3.152 

4.396 

1 

5.805 1 

2.602 

t 

4.258 

5.310 

7.340 

9.651 

D 

1.256 

2.134 

j 2.680 

3.726 

4.910 ' 

1.972 

3.187 

3.957 

5.437 

7.128 

m 

1.152 

1.961 

1 2.463 

3.421 

4.507 ; 

] 

1.698 

2.742 

3.400 

4.666 

6.112 

■ 

1.087 

1.860 

i 2.336 

i 

3.243 

4.273 ! 

1.540 

2.494 

3.091 

4.242 

5. 556 

■1 

1.043 

1.791 

j 2 . 2.50 

3.126 

4.118 i 

1,435 

2 333 

2.894 

3 972 

5.201 

■1 

1.010 

1.740 

j 2.190 

3,042 

4.008 ■ 

1.360 

2 219 

2.755 

3.783 

4.955 

■1 

0.984 

1.702 

I 2.141 

2.977 

3.924 

1.302 

2.133 

2.649 

3.641 

4.772 

10 

0.964 

1.671 

1 2.103 

! 

I 2.927 

3.858 ; 

1.257 

2.065 

2.568 

3.532 

4.629 

11 

0.947 

1 .646 

i 

i 2.073 

2.885 

3.804 1 

1.219 

2.012 

2.. 503 

3.444 

4.515 

12 

0.933 

1.624 

i 2.048 

2.851 

3.760 1 

1.188 

1.966 

2.448 

3.371 

4.420 

13 

0.919 

1.606 

' 2.026 

2.822 

3.722 

1.162 

1.928 

2.403 

3.310 

4.:141 

14 

0.909 

1.591 

! 2.007 

2.796 

3.690 i 

1.139 

1.895 

2.363 

3.257 

4.274 

15 

0.899 i 

j 

1.577 

1.991 

2.776 

3.661 '! 

li 

1.119 

1 866 

2.329 

3.212 

4,215 

16 

1 

0.891 

1 .566 

i 

; 1.977 

2.756 

t 

3.637 ; 

1.101 

1.342 

2.299 

3.172 

4.164 

17 

0.883 i 

1.554 

: 1.964 

' 2.739 i 

3.615 ij 

1.085 

1.820 

2.272 ! 

3.136 ; 

4.118 

18 

0.876 i 

1 .,544 

: 1.951 

■ 2.723 

3.. 595 

1.071 i 

1 1.800 

2.249 ' 

3.106 

4.078 

19 

0.870 j 

1.536 

' 1.942 

1 2.710 ! 

3.577 

1.058 

; 1.781 i 

2.228 ! 

3.078 ; 

4.041 

20 

0.865 i 

i 

1.528 

j 1,933 

1 2.697 

3.561 

,i 

1.046 

[ 1.765 i 

1 ; 

2.208 1 

1 

3.052 ' 

4.009 

21 

0-359 1 

1.520 

1.923 

1 

2.686 

3.545 1 

1 

1 . 0.35 

1 

1.750 ; 

2.190 

1 

3.028 ! 

3.979 

27 

0.854 i 

1.514 

1 1.916 

! 2.675 

3.532 { 

1.025 

' 1.736 ! 

2.174 I 

3.007 1 

3 . 9.52 

23 

0.849 , 

1.508 

’ 1 .907 

i 2.665 

3.520 i 

i.oie 

1.724 1 

2.159 1 

2.987 ; 

3.927 

24 

0.845 : 

1.502 

i 1.901 

1 2.656 ; 

3.. 509 1 

1.007 

! 1.712 1 

2.145 

2.%9 i 

3.904 

25 

0.842 

j 

1.496 

1 1.895 

i 

2.647 

3.497 1 

0.999 i 

1.702 i 

i 

2.132 

2.952 j 

1 

3.882 

30 

0 825 1 

1.475 

1 1.869 

: 2.613 ' 

3.454 i| 

0.966 

1.657 

2.080 1 

2.884 

3.794 

35 

0.812 ; 

1 .458 

, 1.849 i 

1 2.588 

3.421 ; 

0.942 

1.623 

2.041 ; 

2.833 

3.730 

40 

0.803 i 

1 .445 

' 1.834 : 

2.568 : 

3.395 i 

0 923 

1.598 i 

2.010 ' 

2.793 

3.679 

45 

0.795 ; 

1.435 

j 1.821 1 

! 2 . 5.52 

3.375 ! 

0.908 i 

1.577 1 

1.986 1 

2.762 ; 

3.638 

50 

0.788 1 

1.426 

1 i .811 1 

! 2.538 i 

3.358 

0.894 : 

1.560 ; 

j 

1.965 

2.735 1 

i 

i 

3.604 


AdsptMl hy pormiMion from Industrial QunUtif Conlrnl, \ ol, XfV. Mo. 10. April IDTiS, from Jirtirl* mititlod "Tiiblro for On#*-Sidrd 

SluuAticai Tolffnnc* Limits” by (J. j LiCtK*rmsn. 


■*Natrella. M. G., Experimental Statistics . NBS Handbook 91, 1963 
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TABLE 2. lO(Co«tlnu«d). FACTORS FOR ONE-SIDED TOLERANCE LIMITS FOR 

NORMAL DISTRIBUTIONS* 

Tl. iwo nurrca vaiu., l.av. bem corrected to tl» valo« given by U. U. Owen In "factors Trchnt^^ 

. • Plan • Sandlu CorDratlon Monograph SCU-M7, available from llio Clearing nousetor icderal acicnunc and Tcclmical Information. 

U.'pnr-ment'ol Cominctw. Springfield. Va. 22151, The Owen Tables Indicate other errors In the tabic beloiv. not exceeding i m >.io last digit. 


0.75 

0.90 i 

j 

0.95 ' 0.99 

0.999 

0.75 

1 

1 ! 

1 


1 ! 


0.95 0.99 


0.999 


3.804 i 6.158 7.655 ilO.552 13.857 

2.619 4.163 5.143 j 7.042 9.215 


2.149 3.407 4.202 5.741 


9.215 — 

7.501 — 


1.895 

1.732 

1.617 

1.532 

1.465 


3.006 3.707 

2.755 3.399 


5.062 6.612 

4.641 ! 6.061 


2.582 

2.454 


3.188 4.353 


3.031 


2.355 2.911 


1.411 2.275 
1.366 2.210 
1.329 2.155 
1.296 2.108 
1.268 2.068 


1.242 

1.220 

1.200 

1.183 

1.167 


2.032 

2.001 

1.974 

1.949 

1.926 


21 1.152 1.905 

22 1.138 1.887 

23 1.126 1.869 

24 1.114 1.853 

25 1.103 1.838 

30 1.059 1.778 

35 1.025 1.732 

40 0.999 1.697 

45 0.978 1.669 

50 0.961 ; 1.646 


2.815 I 

2.736 

2.670 

2.614 

2.566 

2.523 

2.486 

2.453 

2.423 

2.396 

2.371 

2.350 

2.329 

2.309 

2.292 

i 2.220 
j 2.166 

t 2.126 
j 2.092 
: 2.065 


4.143 

3.981 

3.852 I 
3.747 I 
3 659 
3.585 
3.520 

3.463 
3.415 I 
3.370 I 
3.331 I 
3.295 ' 

3.262 
! 3.233 
j 3.206 
I 3.181 
i 3.158 

I 

I 

: 3.064 
' 2.994 
! 2.941 
! 2.897 
' 2.863 


5.686 

5.414 


5.036 |i 
4.900 i! 
4.787 || 
4.690 1 
4.607 I 

I 

4.534 I 
4.471 I 
4.415 
4.364 'I 
4.319 i 

|l 

4.276 :i 
4.238 
4.204 I 
4 . 171 ij 
, 4.143 :i 


4.022 ! 
3.934 i 
3.866 i 
3.811 i| 
3.766 i 


2.849 

2.490 

2.252 


4.408 

3.856 

3.496 


5.409 

4.730 

4.287 


7.334 

6.411 

5.811 


2.085 3.242 

3.971 

5.389 1 

1.954 3.048 

3.739 

5.075 i 

1.854 2.897 

3.557 

4.828 

1.771 2.773 

3.410 

4.633 

1.702 2.677 1 

3.290 

4.472 

1.645 1 2.592 ^ 

3.189 

4.336 

1.5% 1 2.521 , 

j , 

3.102 

4.224 

1.553 j ‘2.458 I 

3.028 

4.124 

1.514 I 2.405 ■ 

2.%2 

4.038 

1.481 ; 2.357 ! 

2.906 

3.%l 

1.450 ; 2.315 i 

2.855 

3.893 

1.424 1 2.275 ^ 

1 

2.807 

i 3.832 

1 

1.397 ; 2.241 | 

2.768 

j 3.776 

1.376 ; 2.208 i 

2.729 

; 3.727 

1.355 i 2.179 i 

2.693 

i 3.680 

1.336 1 2.154 ' 

2.663 

j 3.638 

1.319 ' 2.129 i 

i 1 

2.632 

' 3.601 

I 

1.249 2.029 1 

2.516 

! 

i 3.446 

1.195 1.957 , 

2.431 

: 3.334 

I 

1.154 1 1.902 

2.365 

j 3.250 

1.122 ! 1.857 

■ 2.313 

1 3.181 

1.096 : 1 .821 

! 2.269’ 

•; 3.124 


9.550 

8.348 

7.566 

7.014 

6.603 

6.284 

6.032 

5.826 

5.651 

5.507 

5.374 

5.26S 

5.16'! 

5.07S 

5.00J 

4.931 

4.86( 

4.80( 

4.75; 

4.701 


*Natrella, M. G. , Experimental Statistics , NBS Handbook 91, 1963 





2-n The Scatter Factor Based on the Tolerance Tnfgrw;>i 

The "scatter factor" is a term which has been used to describe the varia- 
bility in cycles to failure data for purposes of predicting a safe life. It is 
.nalagous to a safety factor. For exao,pIe. if N is the «diao cycles to failore 
at a given stress level, a “design value" N^-N/c where c is a scatter factor 
associated with an appropriate tail probability a. 

As an example, consider the problem described previously, the data of v,ftich is 

given in Table 2.9. The design value is defined as the point above which at 

least P-99% (u-U) of the values are expected to fall with confidence of y-955. 

That value has been established as N^.7.23 thou-cycles. The median cycles to 
failure is 


N » 10^ » 10^'^^ 


18.6 thou-cycles 


(2.39) 


and thus the scatter factor is 


; = N/N, = 18.6/7.23 
=2.58 


(2.40) 


The general form of the 
the expressions of N and N 


scatter factor can be derived by substituting 
in terms of the statistics for Y. Thus 


^ = N/N^ = 10^/[10^^''^a,Y^^] 


(2.41) 


^se of the Preoiction Interval i n Establishing a Design 
Consiaer a random variable Y which is normally distributed ..;ith 
mean and standard deviation u and a (as the above example in Sec 2.10 and 
2.11). A random sample of size n is taken, and the sample mean 7 and sample 
standard deviation ^s .are computed. 
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It is required to make probability predictions about future observations. 
Following is a derivation of the prediction interval and a discussion of how 
it can be used to establish design values as well as provide statistics for 

reliability analysis (see References 14 and 15). 

Consider the distribution of (Y^ - Y) where Y^ is a random variable 
denoting a future observation and Y is the sample mean. Y^ and Y are normally 
distributed and have mean and standard deviation of (u,a) and (y,a/vn) res- 
pectively. Furthermore Y^ and T are independent in a probability sense because 
Y is derived from the first n observations and Y^ corresponds to a future 


observation. Hence (Y - Y) is noma! ly distributed with mean zero and standard 

.... ° O' = o/TT (2.42) 

deviation ° n 

Note that (n-1 )s^/o^ has a distribution with (n-1) degrees of freedom. 


Thus it follows that 




■./TTi 


has a student's t distribution with (n-1) degrees of freedom. (14) 

A rational method of establishing a design value can be based on a prediction 
interval. Let Y^ denote the value below which we expect the next observation 
to fall with probability x. Then 


Yp = Y - G(a,n)s 


(2.43) 


where 

5(a,n) > .'l ♦ I 

where t , is the Students t variate with n - 1 degrees of freedom at proba- 
a,n-l 

bility level a 

*This is a simplified version of the Equivalent Prediction Interval concept 
developed by Wirsching and Hsieh (15): See also Section 4,7. 
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The function G. is shown in Fig. 2,12 for 3 = .01 and .05. As n->-’» the t 
variate approaches the standard normal 2 variate (also shown on Fig. 2.12); 
the radial term approaches one. Thus G approaches z as shown. 

As an example, consider the data of Table 2.9; T= 1.27 and s = 0.127. 
Find the lower 1" prediction interval. For n = 22, G (.01, 22) = 2.57. 

Thus 

Yq = 1.27 - 2.57 (.127) 

= 0.943 

The design value of cycles to failure N is given as, 

Njj • 10°*®^^ - 8.77 (10^ cycles) 


The prediction interval also provides a mechanism for specifying the mean 
and standard deviation for a reliability approach. Use of 7 and s directly in 
a reliability format may lead to significant errors for small sample sizes be- 
cause Y and s are only estimates of the population y and cr. The uncertainity 
in the estimates can be accounted for by expanding the sample standard deviation 
s. 


An alternative expression for the design value is 

= Y - Z 
U a 0 


where is defined as the 'equivalent standard deviation" 
and 2.44, it follows that 


(2.44) 

Comparing Eq 2.43 


<^0 = g(n,a)s (2.45) 

where 

g(n,a) = t 1 / 1 + ^ / z 
a;n-l n a 

= G(n,a) / 2 ^ 


Values of G(n,ot) can be established from Fig. 2.12., v-hich is then divided b'/ 2 
to '^hf.Tin n'n,.'x) 
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SAMPLE SIZE, n 


The value can be interpreted as the expanded standard deviation to 
account for the fact that the estimates T and s are random variables. For use 
in reliability analysis one would state that Y has normal distribution with a 

mean and standard deviation (Y, 

For the data of Table 2.9 and Figure 2.12 

g(.01,22) = G(.01,22)/z 

= 2.57/2.33 = 1.10 
Thus, <j^ * 1.10 (0.127) = 0.140 

Thus we could state in a reliability format that Y-N(1.27, 0.140). 

It is necessary in this approach to specify the "reference value" of a. As a 
general recommendation, the value of a - 0.01 should provide reasonable results. 

2.13 How to Establish a Design Value: Summary Comn^nts 

Given a random sample of data, it is often required to specify a design 
value on the safe side of that data. If the term is a "strength" varaible 
(i.e. design becomes dangerous if variable gets too small), then the design 
value should be on the left side. On the other hand, if the term is a "stress" 
variable (i.e. design becomes dangerous if the variable gets too large) then 
the design value should be on the right side. Thus Yg = Y + Ks. 

There are a number of schemes currently used in industry to specify design 
values for a "strength" variable. These methods are illustrated in Figure 2.13 
using fatigue data presented in Table 2.11. A sumtiary of the methods and com- 
puted values are presented in Table 2.12. 

It is interesting to compare the design values with the actual population 
from which the data is sampled. This particular sample of size n = 10 had 
values larger than average. Thus the lower 1% value (No. 5), based on the 
as-jmption that Uy = Y and ay = s and using normal probabilities, was too high. 
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A SUMMARY OF METHODS TO LOCATE DESIGN VALUE ON 
SAFE SIDE (LOWER) OF DISTRIBUTION 



ACTUAL POPULATION 




Table 2.11 


Data for Example to Illustrate How Design Values 
are Established 


This data on Y was actually sampled from a normal distribution having mean 
Wy = 1.25 and standard deviation ay = 0.12. 


Y. * lognoN. Cycles-to Failure, N. 
’ ^ (103 Cycles) ^ 


1.298 

19.86 

1.351 

22.46 

1.328 

21.30 

1.312 

20.54 

1.095 

12.46 

1.415 

25.98 

1.364 

23.14 

1.232 

17.07 

1.112 

12.95 

1.292 

19.58 


Sample Mean, Y = 1.280 
Sample Std. Dev., s = .1047 


Table 2.12 


An Example of the Various Methods to Establish Design Values 
(Basic Assumption: Y is Normally Distributed) 




The message he^re is that the statistical distribution of Y and s should be 
accounted for in establishing the design value. The prediction interval 
(No. 6) and tolerance interval (No. 7) are the only methods which do this. 

In surima ry it is recommended that the prediction interval or tolerance 
interval be used to define design values. The prediction interval can easily 
be used to provide the data needed for reliability analysis. Tolerance 
intervals are commonly used, but they tend to produce conservative values as 
suggested in Fig. 2.13. 
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Chapte. 3 ANALYSIS OF STRAIN-LIFE DATA: THE LEAST SQUARES LINE 
3.1 Preliminary Considerations 

Statistical analyses of data from a single random variable for the 
purpose of establishing a design value was described in Chapter 2. Examples 
of analyses of cycles to failure data at single stress levels were presented. 
But generally, it 1s required to establish N as a function of S. Therefore, 
tests are performed at various stress levels. 

Typical fatigue data might consist of inelastic strain range versus N 
cycles to failure as illustrated in Figure 3.1. It is necessary to 
characterize this data for design purposes. Two basic methods are employed: 

a. The common approach is to define a design curve on the safe (lower) 
side of the data. 

b. A statistical summary can be presented for a reliability analysis 
or probabilistic design approach. 

A key to strain- life data analysis is the determination of the least 
squares line, an estimate of the median value of N for a given strain (or 
stress) level. This chapter discusses the least squares line and some of 
its characteristics. Chapters 5 and 6 then show how to establish design 
values from strain- life data. 

Least Squares Line (Median Curve Through Strain-Life Pat?) 

3.2.1 Introduction 

Methods for analyzing strain-life Ae-N (or stress-life, S-N) data are 
discussed in this chapter and in the next two. The goal of such data 
analysis is to —'vide a characterization of the Ae-N relationship in a 

form suitable for design purposes. General references on the least 
squares method include References 1, 2 and 3. 
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3.2.2 Median Life Curve 

Consider a constant amplitude strain controlled fatigue test in which 
pairs of data N^.) i = l, n are collected. are the cycles to failure 

associated with strain range and n is the sample size. i--. is the 
independent (or controlled) variable and N is the dependent variable. Data 
from a hypothetical test are shown in Fig. 3.1 plotted on log-log paper. 

Methods of basic linear model analysis are typically used to analy.:e 
fatigue data, and these methods will be used herein to describe life 
relationships. 

Consider first a log transformation of variables. Let 

Y»log N X*log(Ae). (3.2) 

Thus X is the independent variable, Y is the dependent variable. 

Clearly there is no functional relationship between Y and X, but there 
does seem to be some kind of relation. It will be assumed that the data 
(e.g. Fig. 3.1) is a random sample from the following model 

Y(x)*Yq(x) +5 (3.3) 

in which 6 is a normally distributed random variable with mean equal to zero 
and standard deviation equal to a and where 

Y^(x) = a + bx (3-4) 

where a and b are constants. Thus for a specified X, Y is a random variable 
normally distributed having mean and standard deviation 

E(Y| X) = Yq = a + bx (3.5) 

a(Y| X) = a (3.6) 

Note the assumption that a is a constant, not a function of X. The line 
Y^ » a + bx, being the mean of Y, will pass through 'the "center" of the data. 
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Because Y is normally distributed, N :g-|yen will be Icgrcrr^al. 


Thus the median of M, denoted as 1^, is given by Yg^log N. 
original coordinates, the Y line can be written as 

Ae = CN’ 

in which it follows from the above definitions that 


In terms of the 


a = -7 log C 


1/S 


(3.7a) 


(3.7b) 


3.2.3 The Least Squares Line 

The perameters a and b (and thus C andt) and a are not known In advance 
and must be estimated from the data (Se^, n,). i.1,„. £q«t1on 3.2 is used 
to translate the data into (X.. Y.) i-l,n. 

Using the method of least squares, a, b, and a are estimated by a‘, 6, 
and s respectively (3), 

a—, n 


x -)2 


(3.8) 

(3.9) 

(3.10) 


b».i:{x - X)(Y - y-)/ z (X, 

1*1 ^ 1 i«l 1 

a * Y - bf 

Where if and fare the sample means of X and Y respectively. Because each Y. 

IS a random variable, the estimates a. b. and s are also random variables. ' 
The "best fit" line 

Y * a + bX 

(3.11) 

IS called the least squares line. Y is the estimate of Y^ , I.e. the mean of Y 
given x. 

As an example, the PP strainrange data for AIS! 316 presented by 
Saltsman and Halford (4), and shown ir. Fig. 3.2, will be analyzed. This 
data IS presented in Table 3.1, along with statistics associated with this 

no cJeir'""' Povorscd by connrossivo plasticity; 
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Figure 3.2 

The Least Squares Line for Fatigue Data 
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Table 3.1 


Statistical Analysis of AISI 316 PP Data 
[after Saltsman and Halford (4)] 
Sample Size n=7 


Strain Range 
(4e)i 

Cycles to Failure 
Ni 

X^*log(Ae)^ 

Y.*log N. 

.00424 

1700 

-2.373 

3.230 

.00105 

35600 

-2.979 

4.551 

.03508 

120 

-1.455 

2.079 

.03496 

68 

-1.564 

1.832 

.00466 

2333 

-2.332 

3.368 

.02066 

116 

-1.685 

2.064 

.02360 

146 

-1.627 

2.164 


Sample Mean of X 
X = -1.986 

Estimate of a 

a = -.6530 

Other Statistics 
s = 0.1427 

y = 29.663 

'x • six? - . .2934 

* 0.5416 


Sample Meai^ of Y 
Y = 2.746 

Estimate of b 

b = -1.711 
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Equations 3.8, 3.9 and 3.10 can be used to obtain a, b, and s, the 
estimates of a, b, and a, respecti vely. 

a = -0.6530 b = -1.7108 s = 0.1427 (3.12) 

Least squares estimators 8 and C are obtained by inverting Equation 3.7b and 


solving for c and C. 

S . l/b > 0.5845 C = ■ .4152 

A 

Least squares estimators a and b are obtained from Equation 3.7b 


(3.13) 


r * 1/b = 0.5845 

(3.13) 

C = = .4152 

The least squares line, Ae * CN^ is plotted on Fig. 3.2. N is defined as the 
estimate of the median N. 

Note that (a) Y given X is normal and (b) the least squares line is the 
estimate of the mean of yIx.* Therefore, it follows that (a) N|Ae is lognormal 
and (b) the least squares line, N, is the estimate of the median of N|Ae. 

3.2.4 An Alternate Form 

The form of the Ae-N relationship given above is most commonly used in 
SRP literature. However, for probabilistic design purposes, it is more 

convenient to express the "Y^ line" (Equation 3.7) as 

N*k(ae)"’ (3.14) 

Comparing Equations 3.7 and 3.14 the constants K and m in terms of C and t 


are 

K = m = 1/C 

The transformation to linear form involves the substitutions, 
Y » log Ae X »log K a = log K b * m 
The least squares estimates of K and m are 

A. _ A. A 

K =10^ m = b 


(3.15) 
(3*. 16) 
(3.17) 


*Y|X means the random variable Y given a specific value of the random variable X. 
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For the example of Fig. 3.2 

K = .2223 m = -1 .7108 

The estimate of the median curve is 

- 1 7 1 

N = 0.222(As) ' 


(3.18) 

(3.19) 


3.3 Statistical Distribution of the Least Sauargs Estimators^ ^ ^ b 
The estimators a and b will have a bivariate normal distribution. 

expected values are, 

EU) - a = b 

and the covariance matrix is 


The 


(3.20) 


C(a,b) 


(a*/n) 


1 + (^/s^)^ 

-r./s2) 


where X is the sample mean of X and 





2 



-2 

X 


( 3 . 21 ) 


( 3 . 22 ) 


EXAMPLE: Consider the PP strainrange data of Saltsman and Halford (4). 


Substituting appropriate values from Table 3.1 into Equation 3.21, including 
the assumption that s^=a^it follows that an estimate of C(a.b) is 


C(a,b) 


0.042 0.020 

0.020 0.010_ 


(3.23) 


The variance of a and b are respectively 0.042 and 0.010 (standard deviations 
0.205 and 0.100). The covariance of a and b of 0.020 indicates that a and b 
are correlated. The correlation coefficient 
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p 


(3.24) 


_ Cov(a >b) 

a^a/. 
a b 


0.020 

“ (. 205 )( 0 . 10 ) 


0.976 


( 3 . 25 ) 


There exists a high degree or correlation between a and b. This 
dependency leads to operational difficulties when we attempt to present 
fatigue data summaries in a format useful for designers. Practical 
representations of these design factors will be the subject of Chapter 4. 


3.4 Confidence Intervals for the Parameters and the Least Squares Line 
Because a and b are estimates of a and b, they are random variables. 

For example, suppose that a sample of size n*7 of AISI 316 PP data were 
collected in addition to that of Fig. 3.2. Because it is a random sample, 
a different value of a and 6 and a different least squares line will be 

A A 

obtained. Indeed each time a sample of size n is chosen, a different a, b 

and therefore Y will result. This point is important because the uncer- 
0 

tainty in a and b has a significant effect on risk estimates. 


The consequences of having a and b as correlated random variables is 
reflected in the character of the confidence intervals of the least squares 
line. Summarized in this section is how to compute confidence intervals 
for K and m as well as for Vq(x). 

The AISI 316 PP data given in Fig. 3.2 and Table 3.1 is used as an 
example. Table 3.1 summarizes basic statistical analyses. 

The 1-0 confidence limits for a and b are (3) 

*V2-.n-2 


(For a) 
(For b) 
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where ^^/2*n-2 percentage point of Student's distribution 

with n-2 degrees of freedom and using the data of Table 3.1, the 99S 
confidence limits {a = .01) are computed as follows 


’.005,5 


(For a) 


■.6530 + (4.032)(.1427) |l + (-1.986)^ 

sl7 {7)(.2934} 


= -.6530 + .8265 

(-1 .479, +0.1735) (3.28) 

b) -1.711 + (4. 032)(.1427)/£aT. 5416)3 

= -1.711+ .401 

(-2.111 , -1 .310) (3.29) 

The 99% confidence limits can be described by the probability statement, 

P(-1.479 < a < 0.1735) = .99 ( 3 . 30 ) 

P(-2.111 < a < -1.310) = .99 (3.31) 

Because K and m are monotonic functions of a and b, the confidence limits 

can be evaluated by direct substitution. Results are given in Table 3.2. 

Table 3.2 



Least Squares Estimate 

99% Confidence Limits 
Lower Upper 

m 

-1.711 

-2.111 

-1.310 

k 

0.223 

.0331 

1.491 


Because much of the work in SRP has used the form of Ae=C(Ae)*^ it is 
of interest to consider the confidence limits of C and 4 for the AISI 316 
PP data. Results are sunmarized in Table 3.3. 


Confidence limits on Y^, the mean of Y given X, are computed using 

\ *. ‘a 2;n-2 ^ jl . ijSrpf 

ns'^ 

X 


(3.32) 
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For the AISI 316 PP strain data of Fig. 3.2 and Table 3.1, the 99% con- 
fidence intervals on Y are plotted in Fig. 3.3. 

0 


Table 3.3 

Confidence Limits for c and c (AISI 316 PP Data) 



Least Squares Estimate 

99% Confidence Limits 
Lower Upper 

c 

-.585 

-.437 

-.763 



•k 

★ 

c 

.415 

.0743 

1.208 


Values obtained from Equation 3.13 by noting that the correlation 
coefficient between 1 and 6 p*-! . Upper limit on a was combined 
wit3» lower on b and vice versa. 

The purpose of the plot of Fig. 3.3 is to illustrate the fact that when 
dealing with a phenomena that has as much scatter as fatigue, sample sizes 
larger than n=7 are required for a reliable prediction of the least squares 

A A 

line. Assuming that identical statistics (a, b, s, resulted from 

experiments having different sample sizes, an illustration of the effect 
of sample size can be provided. Fig. 3.4 shows 99% confidence intervals 
for Yq for n=7, 15 and 30. 

For design purposes however, it is customary not to use the least 
squares curve to characterize the data, but rather a curve on the safe 
(lower) side of the data. This approach will be pursued in Chapter 4. 
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Figure 3.3 

99 % Confidence Interval for Least Squares Line 
[AISI 316 Data; after Saltsman and Halford (4)] 
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Figure 3,4 

99 % Confidence Intervals for the Least Squares Line 
for N = 7, 15 and 30 Using Same Statistics 
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CHAPTER 4 STRAIN-LIFE DATA ANALYSIS: GENERAL CONSIDERATIONS 

4 . 1 General Considerations for Design 

The general problem of analysis of fatigue data is illustrated by Figure 
4.1. The data typically will be curved even on a log-log plot and will have 
a scatter (standard deviation of N given Ae, denoted as NjAe) which increases 
with smaller Ae. 

A median curve can be established through the data, typically using a 
least squares method. But the median curve only defines the 50% point. A 
low risk design would require definition of a design curve on the safe side 
of the data as shown in Figure 4.1. 

The design curve of Figure 4.1, denoted as N^, is defined as a curve 
below which one would expect the next data point to fall with probability a. 
Given that N is a random variable denoting cycles to failure, can be 
defined by the probability statement 

P[N<NQjAe] = a (4.1) 

for any Ae. This left tail area, a, is the estimate of the fraction of 
points that would be expected to fall below N^. It can also be thought of 
as a "risk level" although it would be equal to the probability of failure 
only in the case where all other design factors, including strain history, 
were exactly known and the factor of safety equal one. In practice, 
the designer must select the value of a. Typical values which have 
been used range from a = 0.1% to 5%. 
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Figure 4.1 


An Illustration of a Design Curve Based on Risk Level a 


STRAIN, A£ 
(or STRESS) 



(log scale) 
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An approach using methods of modern probabilistic design theory does 
not require a design curve. However statistics of the design parameters 
must be given. For example, it will be shown later for a probabilistic 
fatigue design approach that required information from the linear Ae-N 

model. N=K(Ae)f includes the value of m and the median and coefficient of 
variation of K. 

In summary, two approaches to fatigue design can be used. Each requires 

a different way of characterizing Ae-N data but both require statistical 
analyses. 

^ ^^■^’^^n^stic approach in which a design curve (Ae-N) is specified. 

All other design factors are treated as deterministic or constant. 

2. A probabilistic approach in which statistics of the model parameters 

are stated. All design factors are treated as random variables. 

Both methods are described herein. This chapter summarizes the analysis 
for the simple linear model case. The following chapter explores this 
problem for the more difficult non-linear problem illustrated in Figure 4.1 
in which the scatter is a function of strain level. 

4.2 Definition of a Failure Trajectory 

The concept of a failure trajectory has been very useful in analyzing 
fatigue data for design purposes. In Figure 4.2, the star represents a specimen 
which has failed. The dashed lines (failure trajectories) are drawn through 
the point according to some predefined rule. Typically in this case the 
failure trajectories are parallel to the median curve. It is argued that the 
curve defines the cycles to failure of that specimen if it had been tested at 
a different stress level. For example, a specimen which is "weak" at a high 
stress level would also be weak at a low stress level. 
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The assumption that failure trajectories accurately describe material be- 
havior leads to many useful statistical descriptions that otherwise would not 
be possible. For example we can construct the distribution of te given N as 
shown in Figure 4.2. The design curve defines the lower a.% of failures in both 
directions. This statement cannot be made without the failure trajectory assump- 
tion. 

As will be demonstrated in the following sections, not always will this 
concept be employed, but often it can be quite useful. 

4.3 Linear Life Relationships 

It is assumed in this chapter that the strain-life relationship will be 
described by the basic linear model (Equation 3.3) 

Y » a + bx + 5 (4.2) 

where 6 is a random variable having mean value of zero and standard deviation 
a. Let 

Y=log N x=log S a*log K m=b (4.3) 

Define the random variable by the relation 

5 = log K (4.4) 

Combining Equations 4.2, 4.3, and 4.4 yields, 

N * K^K(Acf (4.5) 
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Because <S ^ N(0,a), it follows that is a lognornally distributed 

random variable. The median and coefficient of variation of K are for base 

0 


10 logs. 


^0 = ' 


% 


(4.6) 


(See Appendix 1 for a description of the lognormal distribution.) 

Note that because K is a random variable, N is a random variable. The 

median life, N is Equation 4.5 with the median value of K ; K »1 . Thus 

0 0 


N - K(Ae) 


m 


(4.7) 


In chapter 3 we saw how K and m are estimated by the least squares 
estimators K and iri and that the least squares line 

/\ 

«*K(Ae)"’ (4.8) 

is an estimate of the median curve N. 

For example in the stra inrange partitioning (SRP) method, there are four 
such strain-life curves depending upon the strain type. See Manson, Hal- 
ford, and Hirschberg (1), or Saltsman and Halford (2) for details on SRP. 


"op - 

tensile plasticity, compressive plasticity 

N^p * H(Ae)'^ 

tensile creep, compressive plasticity 

Np^ * K(Ae)P 

tensile plasticity, compressive creep 

l«cc " "(4s)' 

tensile creep, compressive creep 


Least squares analysis in each of the four cases can be performed as 
described above. The median curve for each is simply the relationship with 
the least squares estimators. 
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4.4 Uncertainties in the Linear Model 


Upon employing the basic linear model with the least squares line for 

design, it is important to note that uncertainties stem from four sources (3) 

1. Basic scatter inherent in the material behavior, as described by a 

2 . The sample standard deviation, s, and estimate of a, is a random 
variable. 

/\ 

3. The estimate of a is a, a random variable. 

4. The estimate of b is b, a random variable and is correlated to a. 

As the sample size n becomes large, a-*-a, and and s approaches 

a. For all practical purposes a, b and s could be treated as constants fcr 
large n (typically the assumption would be reasonable for n>50). However, 
because of the expense associated with fatigue testing, sample sizes will 
generally be small. Therefore in general it is necessary to give full con- 
sideration to statistical distributions of these estimators. 

Uncertainties in the estimators significantly complicate the develop- 
ment of a model suitable for design purposes. In the case that the data 
can be described by a linear model, the equivalent prediction interval 
described in this chapter provides an approximate solution. 

This chapter summarizes statistical methods of analysis commonly used 
for linear data and provides commentary regarding the performance of each. 

4.5 Lower 2 ct or 3a Limits 

It is common practice to establish a design curve as a line which is 
parallel to the median Ae-N curve and which lies a distance typically 2 or 
3a to the left of the median curve. The value a is taken to be equal to 
the sample standard deviation, s of Equation 3.10. Figure 4.3 illustrates 
a 2a curve, denoted as N„ . 
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An Illustration of the 2a Criterion 




Use of such a criterion implies a risk level a as follows: 


* 

Distance to the left of the Median 

Implied risk level , a 

2a 

0.0227 

3a 

0.00135 


*Note that a=standard deviation of Y=log N 

This approach to a design curve can produce inconsistent levels of risk. 

Some of the pitfalls are illustrated in the following discussion and example. 

On the basis of the available data, the statistical distribution 
of the next observation Y is described using the concept of prediction 
interval. The distribution of YjX is Student’s t with n-2 degrees of 
freedom* To compute probabilities of Y for a given X, the 100 ot% lower 
prediction interval can be used (7,8) 


a = P[Y(x) < (a + bx) - t, _2 / IlIL ^ (,X:,y 1 

In Equations 4.10, s accounts for the scatter inherent in the physical 
phenomena, the term with the radical sign accounts for uncertainty in a 
and b, and the use of the Student's t distribution (rather than the normal) 
accounts for the uncertainty in s. 

As an example, letting a = 0.227, the same risk as implied by 2a curve 
results in the curve as shown in Figure 4.3. This curve given by 


log N = Y = 
^ a 


^l-a;n-2 ^ 



* 

e.g. 


See Reference 3 or 8. 
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is ODtained from Equation 4.10. The curve is interpreted as follows. 

On the basis of the available data, there is 2.27% chance that the next 
point (Y,X) selected would fall below this curve. This illustration demon- 
strates that (a) the distribution of Y|X is a function of X, (b) because 

lies below the use of the 2a criterion results in risks which are 
higher than advertised. 

The 3o criterion is illustratejJ in Figure 4.4. The curve N is for 

a 

a*0. 00135. Again we see a significant difference between the N, curve 

3a 

which implies a risk of a, and the curve which actually defines a. 

However, as n -► «, a -► a, b b, s -*• a, the estimators being con- 
sistent. Moreover, and the radical term of Equation 4.11 

1. Therefore, for large n, approaches In summary, the dif-^or- 

ence in and of Figure 4.3 (and and N3^ of Figure 4.4) is due to 
the uncertainty in the estimators. 

Another way of looking at the problem is to consider the actual risk 
level a as a function of sample size n for the 2a and 3o curves. The 
curves shown in Figure 4.5 illustrate that the actual risk level depends 
upon the sample size and can substantially differ from the implied risk 
level . 

In suinnary, a 2a or 3a criterion can produce consistent probability 
levels in Ae-N data only if sample sizes are sufficiently large. As a rule 
of thumb, a large sample would typically be n > 50. For smaller samples, 
due consideration should be given to the fact that a, b and s are random 
variables. Proposed in Sections 4.7 and 4.8 are models for consistent 
description of Ac-N data. 


llustration of the 3o Criterion 
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VERSUS IMPLIED RISK 
USING 2a AND 3a CRITERIA 







4.6 Implied Risk in Fatigue Design Curves of the ASME 
Boiler and Pressure Vessel Code 

Fatigue evaluation methods for pressure vessels are included in Section 
III, Division I and Section VIII, Division 2 of the ASME Boiler and Pressure 
Vessel Code (4,5). Both sections contain fatigue design curves. 

The fatigue design curves were developed by applying factors of 2 on the 
median strain range and 20 on the median life (s)- 

The median curve is the Langer model (/)» having the form 

S = + S (4.12) 

e 

where S is defined as 

S = l/2(Ae)E (4.13) 

E being the elastic modulus. K and are the parameters to be determined 
from the data. can be interpreted as the endurance limit. The design 
curve is actually a lower bound envelope. See Figure 4.6. Past experience 
indicates that these factors satisfactorily account for possible size effects, 
data scatter, surface-finish influences, moderate environmental effects and 
time history effects on fatigue resistance* 

The factor of 20 on life dominates the design curve for cycle lives 
typically less than 20,000 (6) as suggested by Figure 4.6. It is of interest 

to examine the actual risk implied by such a criterion. 

Figure 4.7 illustrates how the design curve, and risk level a are 

related. Note that a depends on the amount of scatter in N, but does net. 

For a narrow distribution of N|Ae, a is small and vice versa. 

The relationship between a and the coefficient of variation of N, 

denoted as C.., can be derived as follows 
N 

★ 

For an interesting commentary on the ASflE fatigue curves, see the article 
"Safety Factors in Fatigue Design; Arbitrary or Rational", by J.T. Fong 
and J.H. Smith in Critical Issues in Materials and Mechanical Engineering, 
ASME, 1981. 
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Is Established in the ASME 



CYCLES TO FAILURE 
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a = P[N . Np] 

= P[N - N/20j 


(4.14) 


Assumiruj that N has a lognormal distribution (see Appendix 1), it follows 
that 


X = ^("1og(20)/a^^^) ^4.15) 

where 

*^109 N “ '^0*^34 1og(l + (4.16) 

and where ^(*) is the standard normal distribution function. 

In Figure 4.8 the risk level a is plotted as a function of C., for values 
of commonly observed in fatigue test data. Smooth specimen low cycle data 
typically have from 0,30 to 0.70, but welded joint data for example can 
have values of Cj^ as high as 1.5. 

In summary, the plot of Figure 4.8 clearly shows that application of the 
"factor of 20 on life" curve leads to inconsistent levels of risk. 

4.7 The Tolerance Interval Used to Establish 
a Design S-N Curve ~ 

Examples above illustrated how certain rules for drawing design curves 
can lead to inconsistent levels of risk between data sets. Two methods are 
presented herein for constructing, design curves which are statistically 
consistent. They are a) the tolerance interval method described in this 
section, and b) the equivalent prediction interval (EPI) described in 
Section 4.8. The EPI method can also provide appropriate statistics for 
reliability analysis. 
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FIG. 4.3 


RISK LEVEL* a 


RISK LEVEL AS A FUNCTION OF 
COEFFICIENT OF VARIATION OF 
LIFE FOR A DESIGN CURVf BASED 
ON 1/20 OF MEDIAN LIFE 



COEFFICIENT OF VARIATION OF CYCLE LIFE> C,^ 
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The tolerance limit is established for a single variable as described 
in Chapter 2. However, it may be possible to use the tolerance limit to 
define a design curve from SN data. It is necessary, however, to assume 
the existence of failure trajectories , IKnes parallel to 
the least squares line are illustrated in Figure 4.9. Thus SN data can be 
translated into an "equivalent" random sample of N for any S. The standard 
deviation of Y = log N is constant, with the mean of Y defined by the least 
squares line. 

The tolerance limit at any stress level will lie to the left of the 

least squares line Y, a distance of K s. Thus the design curve is 

^ ot,Y 

No • ? - (4.17) 

An illustration of the design curve based on a tolerance interval is shown in 
Figure 4.10. 

The PP data of Reference 2 is used to provide an example of a design 

curve based on the tolerance interval. The data is shown in Figure 4.11. 

Assume that the decision has been made to define the design curve as the line 

below which no more than a=l% of the population is expected to fall with 

confidence y~9S%\ n=*7. From Table 2.10 K =4.64 The design curve is 

ot,Y ^ 

Y - K s = Y - 4.64(0.1427) 

a,Y 

A. 

» Y - 0.662 (4.18) 

The scatter factor is 

; = lo'^.Y^ , ]g0.662 

=4.59 (4.19) 

In summary, the concept of the tolerance interval can be used to establish 
rational and consistent design SN curves. It is required to assume that failure 
trajectories define material behavior. Moreover, the tolerance interval does 
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Figure 4,9 


The Use of Failure Trajectories to Establish an Equivalent Sairole 
of N at any Stress Level 


Strain Range 
or 

Stress Range 



Cycles to Failure, N 
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Figure 4.10 


Design Curve for 5N Data 
Based on Tolerance Interval 
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Use of the Tolerance Interval to Establish Design Curve [after Saltsman and Halford (2 )] 
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not directly provide the information for a probabilistic design approach. 

These restrictions do not apply to the use of the equivalent prediction in- 
terval described in the next section. 

4.3 A Consistent Method for Characterizing Linear Fatigue 
Data - The Equivalent Prediction Interval (EPI) 

The goal of the analysis described in this section is to develop a procedure 
of fatigue data analysis having the following properties: (a) scatter in data is 
represented, (b) consideration given to the statistical distributions of the 
estimators, (c) probability estimates from the model are consistent and accurate. 

(d) appropriate statistics for a probabilistic design analysis are provided and 

(e) the model is easy to use. The equivalent prediction interval proposed herein 

satisfies these requirements. 

Assume that the strain-life relationship is of the form 

N - Wie)" (4.20) 

For a given :le, K and m (and therefore N) are random variables by virtue of the 
scatter in the data. The value of K in Equation 4.20 is now equivalent to K^K 
in Equation 4.5. 

In a design procedure, two strategies exist. 

1. A deterministic approach in which a "design curve" is established on the 
safe side of the data. Assume that YjX is normal with standard 

deviation of s. For example, to establish a line below which only 1 of 
the points would be expected to fall, let a=.01 and z^=2.33 from normal 
tables. The "design curve" is given as 

log Np = Y - z^_^s (4.21) 
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Nq 1s shown in Figure 4.12 (see also Figure 4. ). Use of this as a 
design curve is subject to the pitfalls as described in Section 4.5 
and Section 4.6. 

But this approach insures a lower 1% bound only in the limit as 
n-«“. For small n, typical of fatigue tests, this method is uncon- 
servative as described below. 

2. A probabi li Stic approach in which a and b (and therefore m and K) are 
treated as random variables. Standc.d probabilistic design methods are 
employed to insure structural integrity. But a and b are bivariate 
normal (3), and therefore m and K appear in the failure function as 
correlated random variables. The problem becoires complicated. 

The distribution of a future observation, defined by the prediction 
Interval is described in Section 4.5. The lower 1% prediction interval is 
also shown in Figure 4.12. Again, we see that use of the curve is mis- 
leading in terms of the implied risk. 

It would be difficult to incorporate the notion of a curved prediction 
interval into a design algorithm. However, it is observed that the equal 
p obability N curve is relatively flat. This suggests that prediction 

\X 

' cervals can be approximated as straight lines, called herein, equivalent 
prediction intervals , and abbreviated EPI. 

Development of the EPI concept is described by Wirsching and Hseih 
(8) and summarized in following discussions. 

Define an equivalent constant standard deviation of Y as 

% = s g(n,a) (4.22) 
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Figure 4.12 

The 1% Njj Curve and the Lower 1?. Prediction Interval [after Saltsman and Halford (2)] 
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where 


g(n,a) * exp[A(a){ln n } 


A(a) 

B(a) = 3.32 - 1 .7a 


(4.23) 

(4.24) 


6 < n < 50; 0.01 < a < 0.15 

Basically the idea is to use the model of Equation 4.2 , but with an expande 
value of s (the value of defined below) so that the curve is shifted 

the left to natch the N curve. ^ ^ 

9 (n.a) > 1 is in essence, an adjustsent factor to s to account tor the fact tha. 

there I's uncerUInty in the estimates of a and b and s. For convenience, a 

Chart of g(n,a) is given in Figure 4.13. 

The model, suggested by the above discussion, is as follows: 

1. Let m » b be a constant. 

2. Assume that all of the uncertainty due to scatter in the data is accounted 
for in a (and therefore K) by considering the y intercept as a random 

variable. 

3. Therefore, let the empirical relationship be 

(4.25) 

Y = a^ + bx 

where a, has a normal distribution with mean ^ and standard deviation o„. 
The consequences of such a model are 

1. Y|x has a normal distribution. (Thus N given Ac has a lognormal distri- 

bution) 

2. The mean value of Y|X is S t bx. (Thus the estimate of the median of N is 

N - ic(Aer 

3 The standard deviation of YlX is (and is not a function of X). 

4 . : . iog K is normal and X is lognormal. The median K end coefficient of 

variation C, of K can be obtained from the lognormal (base 10) forms 


4-25 


Figure 4.13 

FACTOR g(a,n) TO MULTIPLY 
SAMPLE STANDARD DEVIATION 
TO OBTAIN EQUIVALENT 
PREDICTION INTERVAL 




1 


(4.27) 


Ck = _ 

EXAMPLE Given the fatigue data (n = 7) as illustrated in Figure 4,12 

define a design S-N line which is estimated tc be on the safe side 
of 99^0 of the data. This line is to be the a = 0.01 EPI. The basic data 
is summarized in Table 3.1. The estimators are presented in Figure 4.12. 

Using Equation 4.23 with n = 7, a * 0.01 

A(a) = 4,64 B(a) = 3.30 g(n,a) = 1.67 (4.28) 

Thus, the equivalent standard deviation is, 

Oq * g • s (4.29) 

* (1.67)(0.143) * 0.239 

Design Curve The \% EPI is given as N*, where 

log N = Y - ( 4 ^ 30 ) 

This EPI could be used as the design curve in the conventional approach. 
The EPI is shown in Figure 4.14 along with the prediction interval N which the 

OL 

EPI approximates. 

Figure 4.14 suggests that the EPI is a reasonable approximation to the 
prediction interval. As the sample size becomes larger, becomes flatter and 
the EPI becomes an even better approximation [7]. 

Probabilistic Format The data will be analyzed in a format which is convenient 
for probabilistic design procedures. 

The fatigue equation will contain m and K. Using the method described 
above, m is a constant and equal to 

m = b = -1.71 (4.31) 
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(4.32) 


Kwill be lognomal with a median value of (Equation 4.26) 

K = = .222 


Kwill have a coefficient of variation of (Equation 4.33) based on a - 0.01 

- 1 = 0.595 (4.33) 

where 7 was obtained from Equation 4.27, 

0 

In summary. Chapter 4 discussed the problem of establishing a design 
curve and/or appropriate statistics for the case where strain-life data has 
o linear trend (on a log-log basis) and uniform scatter. Analysis associated 
with more complicated strain life models are treated in Chapter 5. 
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Chapter 5 STRAIN-LIFE DATA ANALYSIS; SOME ADVANCED TOPICS 
5.1 Preliminary Remarks 

Chapter 4 discussed how to establishe a design curve, and provide 
parameters for a probabilistic design format for strain-life data which 
a linear trend throughout the range of -'nterest, and has a constant standard 
deviation of log N given log S . This chapter will describe analysis 

procedures for other types of SN data. 

1. The SN data has a linear trend, but scatter increases with decreasing 
stress (or strain) 

2. The general strain life relationship which is now commonly used in design 

(1), i.e., for the local strain approach 

3. A distribution of stress given life N, in the high cycle range where a 

stress endurance limit is assumed 

5.2 Linear SN Curve With Variable Scatter 

Consider the fatigue data of Figure 5.1. Visual examination leads one to 
the conclusion that while there is a linear trend to the data, it does appear 
that the scatter in cycles to failure N increases with decreasing strain e 
(or stress). Data having this characteristic is said to be "heteroscedastic". 
This behavior is fairly common in fatigue data. 

The following analysis procedure is proposed to characterize the data for 
design purposes. The data of Figure 5.1 is used as an example. 

1. It is assumed that (a) the data follows a linear relationship, (b) the 
distribution of Y = log N given X = log Ae is normal. The least squares 

/V /s. 

line Y = a + bx, is obtained from the data (X.,Y.), i = l,n. It is an 
estimate of the mean of Y for a given X and therefore the median of N 
for a given Ae, The least squares line is shown in Figure 5.2. 


5-1 


Figure 5.1 

Strain-Life Data Indicating Increased Scatter at Lower Strain Levels 
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2. Define the deviations as d. 


N 

(Y - Y. ) . Compute, 


s 


1 


n 


> y d . 

n-2i;, 1 


2 


(5.1) 


as an estimate of j , the average standard deviation of Y given X. For 

avg 

the data of Figure 5.1, s = 0.194. 

3. Assume that !dj, the absolute value of the deviation, has a linear rela- 
tionship with X. Determine the least squares line |d(x)| = c + eX. 

This is shown on Figure 5.3. 

4. Assume that the standard deviation of Y given X has a linear relationship 
to X. It is denoted as a(x). The estimate of a(x) is denoted as s(x). 

A 

Assume that s(x) is proportional to |d(x)|. The scale of s(x) is esta- 
blished by defining s (the estimate of a^„_) as the value of s(x) corres- 
ponding to the mean value of X, denoted as T. In the example X » -2.05. 
Both X and s are shown in Figure 5.3. The scale of s(x) is established 
by noting that |d(x)l and s(x) are zero at the same point. 

To use this information for design, several of the techniques described 
previously can be used. As examples, 

1. A design curve defined as the mean minus M standard deviations can be 
easily constructed from Figure 5.3. The M » 3 curve shown in Figure 5.2 
illustrates the heteroscedastic nature of the data 

2. Assume that fatigue behavior can be described by failure trajectories 
which are flared from the least squares line in proportion to |d(x)|. An 
example is shown in Figure 5.4. Then, for example, a lower tolerance 
limit can be used as a design curve as shown. The values of P = 0.01 and 
Y = 0.95 are those of the "A - values" in MIL-HDBK-5. 
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Least Squares Line and Design Curve for Heteroscedastic Data of Figure 5.1 
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3. For a reliability format, it is necessary to have the least squares 
line, and the standard deviation. The standard deviation has the 
functional form, 

a(x) = A + BX (5.2) 

For example, the estimate of a(x) is, from Figure 5.5, 

s(x) * -0.0834 - 0.1377X (5.3) 

This representation of a(x) significantly complicates reliability analysis. 

However, a simplification using an equivalent standard deviation c can be made. 

o 

Assume that a variable amplitude loading, consisting of k discrete strain ampli- 
tude levels, e^, i * l,k are applied. The fraction of time that t. is applied 
is known to be f.. An average, or equivalent standard deviation can be defined 
as a weighted average, ^ 

^0 (5.4) 

1=1 

where X^. = log e^. But the right hand side is, by definition, the expected value 
of c(x), denoted as E[a(x)]. Substituting Equation 5.2, 

Oq = E[a(x)] = E[A + BX] 

Oq = A + B*E[log e] (5.5) 

To find E[log e] it is necessary to evaluate 

k 

E[log e] =_^^f Jog e. (5.6) 
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and if the strain amplitude levels are continuous with probability density 
function f(c), then 

E[log e] = [ (log t)f{z)dt (5.7) 

- 0 

In summary, the value of as computed by Equation 5.5 can be used in a 

Av A 

reliability analysis. The least squares, NS" = K defines the median value of 
K, and the coefficient of variation of K is, (using lognormal relationships) 

C|^ - 1 (5.8) 

5.3 Statistical Analysis of the General Strain-Life Relationships 

Fatigue tests can include cyclic plasticity and accompanying low cycle 
fatigue. A stable hysteresis loop is recorded for each specimen as shown in 
Figure 5.5. The figure also defines strain types. Strain reversals to failure, 
2N, is recorded for each specimen. The general strain life relationship is 
commonly used to describe fatigue strength for the modern local strain approach 
to fatigue. General references include References 2 and 3. 

An example is presented in this discussion. The data as well as the 
analyses are given in Table 5.1. 

The Coffin-Manson rule for the plastic strain-life relationship is 



The empirical constants (fatigue ductility coefficient) and c (fatigue 
ductility exponent) are established by a least squares line through the data 
as illustrated in Figure 5.6. 
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total strain range 
plastic strain range 
ACp/2 = plastic strain amplitude 
: elastic strain range 
: Ac = elastic strain amplitude 




Strain-Life Data for Example Problem 
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Estimates of the standard deviation of log e (vertical deviations) 


PLfl|TIC ST 



(5.10) 


The Basquin equation for the elastic strain-life relationship is, (4) 


~ea ■ 2 

Again, the empirical constants S'^ 
method. See Figure 5.7. 

The total strain life curve is 
strains . 



and b are established by a least squares 


obtained by adding the elastic and plastic 


r- ^ 


+ £ 

ea pa 


’ Y (2N)^ + (2N)'^ (5.11) 

In this example, the total strain-life curve is shown on Figure 5.8. 

The "standard deviations" in Table 5.1 were computed as follows. Consider 

A 

A A 

a failure point (c^, N^. ). Let ^i * ^°^10^i ^i 

value strain from the least squares line at N.. Let 

d. - X. - X. (5.12) 

The sample standard deviation, an estimate of is given as 


s = 


_L y d 

n-2i,i 1 


n 

2 


1/2 


(5.13) 


Values of s, so computed, are given in Table 5.1, not only for the plastic and 
elastic strain curves, but also for the total strain curve. 

For design purposes, it may be required to define a design curve on the 
safe (lower) side of the data and/or specify statistics of the parameters. The 
latter would be necessary if a complete fatigue reliability analysis is to be 
performed. 

At this time only preliminary studies of the statistical character of 
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strain-life data have been completed. Analyses of 14 sets of fatigue data 
of steels and nickel base superalloys of various types have suggested the 
following conclusions. 

1. There is no statistical relationship between the plastic and elastic 
strain amplitude for a failed specimen. One might suspect one of 
two possibilities (a) a specimen having high ductility would have 
low strength,... good plastic strength, poor elastic strength, or 
(b) a defect causes poor strength in both plastic and elastic strain. 
No evidence in these preliminary studies suggests either phenomena 

is dominant. 

2. Scatter in cycle life (N) increases with decreasing strain level (e,). 

a 

This behavior is typical of fatigue data. 

3. Scatter in fatigue strength (c.) decreases with increasing life (N). 
Experimentalists have suggested to the author that tie standard devia- 
tion of log e. is constant (not a function of N). No evidence that 

d 

this is fact is found in the present study. 

The above conclusions are preliminary, and work on this problem is 
continuing. 

Based on the observations as described above, the following procedures 
characterizing the data are presented as being reasonable. 

1. The lower a% EPI curves associated with the plastic and elastic strain 
life data can be established as described in Chapter 4. The 1% EPI 
curves are shown in Finures 5.6 and 5.7. Then the total 1% EPI curve 
of Figure 5,3 is obtained by adding (vertically) the two. 

2. A lower tolerance limit can be established in the same way, i.e. by 

first using conventional methods (see Chapter 4) for e ,-N and e -N 

pa ea 

curves and then adding the two. The result is shown in Figure 5.9. 
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Ffgure 5.9 

Strain-Life Curve and Tolerance Limit for Example Problem 
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3. For a reliability analysis, statistics of Sp b and c are required. 
Assume b and c are constants and equal to the least squares estimators. 
Assume that s^ and are independent random variables having lognormal 
distributions. In the example of Table 5.1, the statistics required 
for a reliability analysis are as follows. 


Table 5.2 

Statistics for Reliability Analyses for Example Problem 



Median 

Coefficient of Variation 

»f 

148.8 ksi 

0.042* 

^f 

.491 

0.263* 

b 

-.108 

★★ 

0 

c 

1 

1 

cn 

o 

** 

0 

based 

on a 5S EPI , 

these values would be only s 


from the IS EPI values given in Figures 5.6, 5,7 and 5.8 

'Ar 

variable not random 
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5.4 A Model for the Distribution of Fatigue 
Strengths at High Cycle Live? 

5.4.1 Preliminary Remarks 

The number of lead cycles experienced during the lifetime of a wide 
variety of structural and mechanical components is in the range from 10^ 
to 10^^ cycles. When metal fatigue is considered to be the mode of failure, 
it is necessary for the designer to know the fatigue strength (the maximum 
amplitude of the oscillatory stress) for a given cycle life fJ. 

Fatigue testing is however expensive, partly because of the length 
of time that it takes to apply these long cycle lives. An accelerated test 
method which has been used successfully involves loading of a specimen at 
a relatively high stress level, with relatively short cycle lives, and 
then extrapolating the data to the higher cycle lives. 

Described as follows is a study of an analytical method for the 
statistical distribution of fatigue strengths at high cycle lives based on 
data at relatively low cycle lives. 

5.4.2 Basic Assumptions and Description of the Problem 

Consider a constant amplitude high cycle fatigue test in which pairs 
of data (s!,N.), i = 1, n are collected. N. (the dependent variable) is 

I 

the cycles to failure associated with stress amplitude (the independent 
variable) and n is the sample size. Data from a hypothetical test are shown 
in Figure 5. ID. An endurance limit is assumed as illustrated by the horizontal 
segment of this S-N curve. It is assumed that the endurance limit occurs at 
10^ cycles. Furthermore, it is assumed that the distribution of S at a given N 
at 10^ cycles will be the same for any higher life cycle as shown in Finure 5.10. 
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A MODEL FOR HIGH CYCLE FATIGUE 
(Typical Use; Machine Components) 
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In general, a fatigue test plan specifies a termination of the 
loading on a component which does not fail before a given life. This spec- 
imen is called a "runner" and the data points are shown by the triangles 
of Figure 5. IQ in which the test was suspended at 10^ cycles. This data 
must be included in the analysis. 

To obtain a random sample of S, now defined as a random variable 
denoting the fatigue strength at 10^ cycles, the following procedure is 
used. It is assumed that the fatigue strength for a given specimen is 
defined by a straight line (failure trajectory) connecting the fatigue 
strength coefficient at s^ 1/2 cycle and (N^,s^) as shown in Figure 5.10. 

The data point (solid point in Figure 5.1?) is projected to the life cycle 
where the endurai.ee limit occurs, which is assumed here to be 10^ cycles. 

The sample point (circle in Figure 5.10) is denoted as S^. By such a 
scheme one can obtain a random sample of s. 

5.4.3 Analysis of Suspended Fatigue Test Data 

Incomplete failure data consisting of levels on failed components 
and unfailed components are called "multiply censored." This suspended 
data can be analyzed utilizing the "median rank concept" and suspended 
items approach [5]. The median ranks extracted from this approach will 
be used to establish the empirical distribution function of fatigue 
strengths . 

Lipson and Sheth have described a method for obtaining the empirical 
distribution function from a random sample which includes suspended items [5]. 
This method is a combination of a modified sudden death approach and Johnson's 
concept of median ranking [e]. The method involves an adjustment of the order 
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number based upon the position of the suspended items. Step-by-step instruc- 
are illustrated by tie followinq; 

1. Suppose that the failure data on n specimens consist of the failure 
stresses (at 10^ cycles for the failed units and the stress levels for the un- 
failed units (see Table 5.3 and Figure 5.11 ). Ordrr the failed units in the 
sample from the smallest to largest failure stresses as shown in Column 3. 

2. Obtain the number of suspended items which precede each failed 
unit (Column 2). 

3. Determining the "new increment" of each failed unit by using the 

formula 

New Increment - 1" * , U : failure order "umber 

1 + nunioer of items following present suspended set 

The new increment is recalculated each time a suspension is encountered in the 
ordered stress table (see Column 4), 

4. Calculate the order number for each failed unit. This calculation 
is done by simple addition of the last order number and the new increment 
(Column 5). 

5. Obtain the median rank (empirical distribution function) of the 
order number for each failed unit by the formula (Column 6). 

The empirical distribution function of Column 6 is used as 

a basis for selecting a suitable statistical model. The method of least 
squares, as described in Sections 2.5, 2.6, and 2.7, can be used. The cor- 
relation parameter for the three distributions considered are 
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Table 5.3 

Analysis of Suspended Fatigue Data 


EXAMPLE : 


Given a sample of fatigue failure stress at 10^ cycles (in ksi); 
n = 13; the data are ordered. 


11.7 i 12, 

i 

2| 12, 5| 

12. 7i 12. 8| 

13. 3| 

13.8(14.0* 

Il4.0*jl4 

t 1 

.0*il4.4 ; 14. 7j 15.0* 

♦Specimen did not fail at 10^ cycles 

• 



Organize data as shown below: 





(1) 

(2) 

(3) 


(4) 

(5) 

(6) 

Failure 

Suspended 

Stress 


New 

Order 

Median Rank 

Number 

Items 

S/.x 

Increment** 

Number** 


Preceeding Failure U; 



(j) 

n + 0.4 

1 

0 

11.7 


1.00 

1.00 

0.052 

2 

0 

12.2 


1.00 

2.00 

0.127 

3 

0 

12.5 


1.00 

3.00 

0.201 

4 

0 

12.7 


1.00 

4.00 

0.276 

5 

0 

12.8 


1.00 

5.00 

0.351 

6 

0 

13.3 


1.00 

6.00 

0.425 

7 

0 

13.8 


1.00 

7.00 

0.500 

8 

3 

14.4 


1.75 

8.75 

0.631 

9 

0 

14.7 


1.75 

10.50 

0.761 


**An example of how the new increment and order number are calculated: 


The eight failure is preceded by three suspended items. Therefore, to 
find the increment, as shown in Step 3 above. 

Thus, the new order number is j = 7.00 + 1 .75 = 8.75 

The remaining order numbers are determined by using the same procedure and 


S-N DIAGRAM FOR EXAMPLE OF DATA EXTRAPOLATED INTO HIGH CYCLE RANGE 
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CYCLES TO FAILURE. 


Distribution Family 

Correlation Parameter, c* 

Normal 

0.2010 

Lognormal 

0.1809 

Wei bull 

0.2738 


The lognormal is considered to have the best fit because it has the smallest s. 
The probability plot is given in Figure 5.12. The line is fit by the least 
squares method and the distribution parameters on the figure are computed from 
the least squares line (using Equation 2.15 and 2.16). 

In summary, for example, the distribution of fatigue strength S given 
cycle life N in the high cycle region (beyond 10^ cycles) can be modelled as a 
lognormal distribution. The parameters are = 2.67 and = 0.107, where 

X = In S. Using the relationships of Appendix 1, the sample mean and stan- 
dard deviation of S is S = 13.8 ksi and s^ = 1.48, respectively. Full con- 
sideration was given to "runners" in the data analysis. 

Statistical methods of estimating the fatigue strength at a given cycle 
life are available. Such methods, suiimarized by Reemsnyder (7), Lipson and 
Sheth (5), and by Collins (8) include the survival method, the staircase method, 
and the Prot method. Because these methods typically require long cycle lives 
and relatively large sample sizes, testing tends to become expensive. 


★ 

See Equation 2.29 
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1 1 1 1 1 1 1 1 s— I 1 1 

2.40 2. 45 2.50 2-55 2.60 2.65 2.70 

X = LN( STRESS) 

Estimated Parameters: 

Mean of S = 13.15 
Standard Deviation of S = 1.01 
= 2.62 
a = 0.107 

X 
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CHAPTER 6 LIFE PREDICTION EQUATIONS USED FOR FATIGUE DESIGN 
6.1 Preliminary Remarks 

Engineering models used by designers for fatigue failure avoidance usually 
employ the concept of "damage". Summarized in this chapter are some of the var- 
ious forms of damage in current use. It is shown that the concept of damage can 
be employed for variable amplitude loading using both a characteristic S-N curve 
approach and under certain assumptions the fracture mechanics approach. More- 
over the concept is used also for strainrange partitioning in which damage is 
identified not according to stress level but with each of four different types 
of inelastic strains. 

The limit state is achieved when damage D attains a critical level A. 
Experimental results in which damage at failure A is measured are summarized 
herein. Both variable amplitude loading and strainrange partitioning results 
are presented as statistical distributions of A. 

6.2 The Concept of Damage 

Assume that the SN curve for a component subjected to constant amplitude 
loading is as shown in Figure 6.1 and that no statistical scatter exists. Blocks 
of constant amplitude loads are now applied as shown. Fatigue damage D associated 
with this variable amplitude process is defined as 



where n., and N(S^) and k are defined on Figure 6.1. 

More generally as a component accumulates a load history in service, damage 
can be expressed as a function of time i and several other design factors denoted 
by the vector U, which relate to the loads as well as to the SN curve. Thus 

DAMAGE = D(n, U) or D(x, U) (6-2) 
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Miner's rule states that failure occurs when D, a monotonical ly increasing 
function of time, equals unity. A more general definition is provided in Section 
6.5. The following two sections summarize some damage models in current use. 

6 . 3 Some Expressions for Damage Used by Designers for Variable Amplitude Loadin g 

Following are examples of expressions for fatigue dama^je under variable am- 
plitude (random) loading. The model to be used depends upon how the load history 
and the SN relationship are specified. 

CASE J Often field measurements result in a variable amplitude stress range his- 
togram as shown In Figure 6.2. The fraction of time that the stress range Is at 
level is f^. Thus n^ can be written as a fraction of the total number of 
applied cycles, n 

n^ » f^n (6.3) 

and the damage at cycle life n is. 



CAi'*^ 2 Assume a linear S-N curve, N » KS'"’. Note that the value of K depends 
upon whether stress amplitude or range is considered. Equation 6.4 becomes 

( 6 . 5 ) 

1*1 

But the summation is recognized as the expected value of s"* for a discrete random 

variable (1). ^ 

E(S"’) * I f.S."’ (6.6) 

i*l ^ ^ 

Thus D = ^E(S"’) (6.7) 
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Note that in the special case of constant amplitude loading E(S ) = S, and 
Equation 6.7 becomes 

0 , n s" (6-3) 

IN 

An equivalent constant amplitude stress can be derived by comparing Equation 6.7 
and Equation 6.3 

Sg = [E(S"’)] (6.9) 

This stress is sometimes called "Miner's stress". Assuming Miner's rule works. 

the value of S can be entered into a constant amplitude S-N curve to predict 
e 

life N under variable amplitude loading. 

CASE 3 Assume that the distribution of stress ranges S can be treated as a con- 
tinuous random variable, the pdf of which is shown in Figure 6.3. The fraction 
of cycles at stress level s in the interval (s, s As) is 


fi » fs(s)AS 


( 6 . 10 ) 


Combining Equation 6.1 and 6,3 with Equation 6.10 to obtain the continuous 
equivalent of Equation 6.4, 


D - n y %(s)As 
i=l N(s) 


( 6 . 11 ) 


In the limit as As 0, 

^ 1 ^c(s)ds 

“ ■ "J -hr- 

0 

CASE 4 When S is continuous and it can be assumed that the S-N curve is linear 
of the form NS^ = K, the expression for damage becomes 
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^ V 


0 = -Hs r5is)as (6.13) 

i^ote that the integral is by definition the expected value of s'”, denoted as E(s’”). 
Thus, as in the discrete case, 

0=fE(s"') (6. 14) 

Any statistical mjdel for S can be used, but it is common to use the Weibull, 
having a distribution function of the form 


Fj(s) « 1 - exp[-(s/5)^] (6.15) 

where the parameters C and 5 must be determined from an analysis of the stress 
data. For the Weibull 

£( 3 ") . 1 ) ( 6 . 16 ) 


In the special case where ^ = 2, the Weibull reduces to the Rayleigh dis- 
tribution. This is an important case because the Rayleigh is the distribution 
of peaks or ranges in a stationary narrow band gaussian process having RMS value 
of a. 


5 = 


(✓7a) if stress amplitudes are considered 
(2^7a) if stress ranges are considered 


Thus we get the familiar form of Miles' eguatlon for damage under n cycles 
of stress of o narrow band process In the case where K is based on amplitude (2). 


D » (n/K) ( i/?a)*''r(m/2 + 1) 


(6.17) 
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6.4 Damage Using a Fracture Mechanics Approach to Fatigue 

Consider the Paris equation for fatigue crack propagation under a sequence 
of constant amplitude stress ranges, S 

da/dn = (6.18) 

where C and m are empirical constants, a is the crack length and n is cycles, 
is the range of stress intensity factor 

A< * YSi^ (6.19) 


where Y is the geometry factor. 

An expression for cycles to failure can be derived by integrating Equation 

6.18, a from a^ (the initial crack size) to a^ (the critical crack length), and n 

from 0 to N (the cycles to failure). Assuming that Y is constant (not a function 

of a) and a > > a^, it follows that 
c 0 


MS'" 


1 

a^^^ ■ ^m/2 - DCY*" 


( 6 . 20 ) 


But this form is identical to the form of the SN curve used in the previous section. 
The right hand side is identified with K, an empirical constant. 

Consider a random load sequence of S^, i * l,n. if sequence effects were not 

important, it may be reasonable to use Miner’s rule as above, and express damage as 

O.Je(s“') (6.21) 


For this to oe a useful inode 1 , there should be a reasonable distribution of S.'s, 
e.g. no grouping of the large and small amplitudes. There should be no large 
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;v.e;'1-aa: to cause significant crack relarcfation. tXor'eove.-', 't "crt be assctSu 

that the mean stress associated with each cycle does not influence the constants, 
C and m. 


As described in Section 6.3, a Miner's stress S can be defined by Eouation 

e 

6.9. Then an equivalent stress intensity factor can be defined as 


= YS-.T (6.22) 

e e 

where = [ECS*") 

If the model is valid, a da/dn < A< curve observed from random loading 
using Equation 6.22 should produce the same results as the constant amplitude 
curve. 


6.5 Damage Exenressions for High Temperature Low Cycle Fatigue 

The TcUajx Lou/: A damage expression which includes the effects of fatigue 
★ 

and creep is 


D * 



(6.23) 


where definitions of n. and are given in Section 6.2 (see Figure 6.1); t^. is 
the time duration of the load condition, i ,and T^. is the time to rupture at load 
condition i. There are q discrete conditions of creep. 

This model is easy for designers to use, but it does not account for syner- 
gistic effects of fatigue and creep. 

StvUnAangz ?(vxtitio*u.ng (SUP} : The method of strainrange partitioning (SPP) 
developed by Manson, Halford, and Hirschberg [3] is well documented in their 
report as well as the later reports of Hirschberg and Halford [4] and Bernstein 
[5]. 

Hr 

A combination of Miner's rule anl p^ohinson's law. 
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Four different types of strainranges can be defined as the basic building 
blocks for any conceivable hysteresis loop; 

1. PP — tensile plasticity reversed by compressive plasticity 

2. CP — tensile creep reversed by compressive plasticity 

3. PC — tensile plasticity reversed by compressive creep 

4. CC — tensile creep reversed by compressive creep. 

First strain-life relationships for each of the strainrange types are 
established by experiment 

N. = A.(Ae)*^' i = PP, CP, PC, CC (6.24) 

Given a hy.;teres1s loop for fatigue stress at a point (obtained from stress 

analysis), the fraction of each strainrange type f^, a component of the total 

inelastic strainrange is identified using an algorithm as described in References 

3 and 4. For example, f ^ where Ae„ is PP strainrange and Ae. 

pp pp 'n pp in 

is the total inelastic strainrange. 


i£. = y Ae. 

1'’ 1=1 


4 


Damage after application of n cycles 
is defined as 

4 

D = n I(f./N.) 
1=1 T T 


i = PP. CP, PC, CC (6.25) 

of the constant amplitude strain cycle, 

i = PP, CP, PC, CC (6.26) 


6.6 Damage at Failure— Some General Remarks 

Assuming that no healing occurs in the component, fatigue damage will be 
a monotonically increasing function of n (or t) as shown in Figure 6.4. Dam- 
age D will increase until it reaches a failure level A. 
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Figure 6.4 

Damage as a Function of Time 



A more advanced model using fracture mechanics methods, considers 
strength as a process which monotonical ly decreases: See "Relia- 
bility Analysis of Fatigue Sensitive Aircraft Structures Under 
Random Loading and Periodic Inspection" by J.rj. Yang and W.J. 
Trapp, AFf^-TR-T^-ZQ , February 1974. 
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Note that if only one stress level is applied, failure occurs when D = I = 1. 
Without any knowledge of the physical process of fatigue, the obvious assumption 
to make is thi-t in general failure occurs when D = A = 1 for any sequence of 
stresses. Indeed this is the Palmgren-Miner hypothesis, or "Miner's Rule" as it 
is often called [6»7]- But recent studies have suggested that for design pur- 
poses it is useful to model A as a random variable [8]. 

Consider a hypothetical fatigue test in which the SN behavior is determin- 
istic as shown in Figure 6.1. Damage at failure denoted as is measured for 
m specimens subjected to different load histories as shown in Figure 6.5. If 
Miner's rule is exact, all = 1. But Miner's rule is a simple model used to 
describe a complicated phenomena. As a result, primarily of sequence effects, 
one should expect statistical scatter in the sample A^, i » l,m. Therefore 
damage at failure A should in general be treated as a random variable. 

In practice of course, observed values of A. will also contain scatter by 
virtue of basic variability in material behavior. To measure the statistics on 
A, somehow material variability effects must be removed. 

6.7 Fatigue Design Relationships 

The design life in cycles (or T^ in time) is specified. Three equivalent 
formulations are commonly used for design purposes. 

fonmilation. J; Damage at the design life N is denoted as D(N^,U). For a 

model in which the designer assumes that A and U are constants, the condition 

for a safe design at life N is 

0 

D(Nq,U) < A (6.27) 
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To insure a safe design given statistical scatter and uncertainties, the 
designer can use a) a design curve on the safe side of SN data, as well as 
"safe-side" values for all U. and/or b) a reduced value of a, e.g., A = 0.30, 

rather than the A = 1 of Miner's rule-. 

Assume that A and U are treated as random variables. The density func- 
tion for each at a given life is shown in Figure 6.4. The probability of 
failure at life n can be written as, 

* P[A < D(Nq,U)] (6.28) 

Example: Consider the damage expression when stress is a stationary 

narrow band process (Equation 6.17). The probability of failure is 

Pf = P[A < (N^/K)(/2a)'"r(m/2 + 1)3 6.29 

where in general A, K, a and m could be random variables. The problem of 
evaluating p^ is described in Chapter 7. 

VoKMiliVtion 2: Let N be a random variable denoting cycles to failure. Thus 

D(M,U) is damage at failure, A 

A = D(N,U) (6.30) 

Inverting the expression, 

N = N(A,U) (6.31) 

Assuming that A and U are constants, the condition for a safe design is 

Nq < N (or T^ < T) (6.32) 
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Given statistical scatter and uncertainty, a safe design can be insured by 
using a) a design curve on the safe side of the SN data, and/or b) conserva- 
tive (high) values of stresses, and/or c) a value of larger than the 

★ 

expected service life For example, = 2N^ has been used [9], 

In a reliability approach, the probability of failure can be written as, 

P. = P[M < N ] (6.33) 

' 0 

= P[N(A,U) < N.3 

-w U 

In which A and U are treated as random variables. Figure 6.6 illustrates the 
pdf of N, as well as P^. 

Example; Consider the damage expression of Equation 6.17 with n = N and 
D = A. Solving for N and substituting into Equation 6.34, 

Pf = P[AK/(v7a)"’r(m/2 + 1 ) < N^] (6.35) 

where in general A, K, o and m would be considered as random variables. 

Formulation 3. Another approach commonly used requires knowing the distribu- 
tion of fatigue strength S at a given design life N . The pdf of SlN is 
illustrated in Figure 6.7. Miner's stress (Equation 6.9) can be compared 
to strength. Miner's stress will be a function of random design factors U. 

In the special case of constant amplitude stress at level S , S = S . 

0 e 0 

If the SN curve is given S * S(N), then the condition for a safe design 
is 

S < S(N ) 
e 0 

applied to an S-N curve already on the safe side of the data 
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For example, if the S-N curve is = K, then a safe design occurs when 
S, < (K/M But, in general, S will be a random variable as suggested 

by the pdf of Figure 6.’. The probability of failure is 

= P[S < Sg(U)] 

Example; Miner's stress for the case of a stationary narrow band Gaussian 
process is obtained by combining Equation 6.9 and 6.16 

Sg » /7[r(m/2 + 1)]^^'" a (6.37) 

Assume that fatigue strength at is given as » K. Thus the probability 
of failure is 

» P[(K/Njj,'^"’ < /?{r(m/2 + 1)}^^"* cj] (6.38) 

In general, K, o and m will be random variables. 

Note that this form is similar to the two above examples with the exception 
that ^ is missing from this form. It's influence could be accounted for by a 
variable which multiplies S^. Such a variable would describe the uncertainty 
associated with using Miner's stress as a characteristic stress. 

6.8 Statistical Considerations of the Palmgren- 
Miner Fatigue Index 

The Palmgren-Mi ner linear damage accumulation rule (referred to herein as 
"Miner's Rule") for estimating variable amplitude fatigue life is easy to use 
and therefore widely employed in design procedures. But because fatigue is 
a complicated process involving many factors. Miner's rule, a simplified 
description of fatigue, does not provide consistently accurate predictions. 
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There is a need to quantify the performance of Miner's rule to provide guid- 
ance to designers. 

A statistical study of the results of fatigue experiments of several 
investigators was performed by Wirsching, Yao, and Stahl, and summarized in 
Reference 3. A suiisnary of this study is presented here. 

Composite statistical models of fatigue damage at failure (denoted 
as d) were developed. The lognormal distribution was used in the statistical 
analysis of A. Such a study of the statistical variability of a relatively 
simple theoretical model is necessary and desirable for the development o>f 
probability-based fatigue design codes. 

Define D as fatigue damage as computed by a linear damage rule. Consider 
a fatigue test performed by an investigator. Let r be the number of specimens 
in the test. Let Ap be the value of D at failure for the pth specimen. To 
evaluate A * 


A = Cycles to failure for the pth specimen 
P Cycles to f^ailure predicted by the PM rule 


(6.39) 


A sample from an investigator would be A^; p = 1, r, with the tests in general 
being performed at different stress levels. 

Because of the scatter observed in the data, it is suggested that the 
quantity A be represented as a random variable. In general, failure can be 
defined as the event (D > A). This is a generalization of Miner's rule which 
states that failure occurs when D * 1. The probability of failure is given by 


= P(D > A) (6.40) 

Assume that A has a lognormal distribution with mean u^, standard deviation 

and a coefficient of variation (C.O.V.) C, = a /u . 

A A A 
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Test data on da;nage at failure i obtained by several investigators Is 
sumr’arized by the distribution functions of Figure 6.8 and in Table 6.1. 
Following is commentary on this data summary and how it was developed. 


The sample mean of A obtained by the ith investigator is 

1 " 


A. 


A. 


and is an estimate of 


' ^ p-1 P 

The sample standard deviation is 


(6.41) 


I ‘"p 




1/2 


( 6 . 42 ) 


and is an estimate of The estimate of the C.O.V. and median of the data 
of the ith investigator are, respectively. 


Ci » s^/A^. 


A^ » A^/ V 1 + C.‘ 


( 6 . 43 ) 

(6.44) 


The assumption that a is lognormal is based on observations of the empi- 
rical distribution functions of Figure 6.8. 


6.8.1 The SAE Fatigue Program 

The SAE Cumulative Fatigue Damage Test Program was an effort to exchange 
information between experts in universities and industry on methods of testing 
and life estimation [18]. The load histories defined were considered to be 
typical of those experienced by components of automotive vehicles. Test 
specimens of RQC-100* and f-lan-Ten** steel were compact tension specimens with 


*Bethlehem Steel; roller quenched and tempered to 100 ksi miniminum 
yield strength. 

S. Steel; hot rolled to 50 ksi minimum yield strength. 
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DISTRIBUTION FUNCTION 






Table 6.1 


1 


Estimacss of Statistics of Jamaqe at failure, ^ 

3ased on a Lonnopmal '•odel for i. 

(n • numoer of specimens) 


.•^an of 0 Median of i 


Standard 
3aviat1on of 




s 


i 


Coef*^icient of 
Variation of \ 


1. Crichlow et a1.(iC );an amalgamation 


of test data of T,any configurations 
and Ttaterlals of flight vehicle ' 1.53 
structures Including full scale ?-51 
and C-46 airplane wing test (n«266) 

Z. SchIJve ( 1 ]J; amalgamation of results 

from 19 reports; all specimens Z024 1.72 

and 707S aluminuai 

3. Richmrt and (leMnarlc CL2); axially 

loaded large (3 1/2 x 1/2 In) 1.27 

specimens of ASTM-A7 (n*31) (K79)^ 

4. Rlchert and M em w ark (12); rotating 

bean small specimens (dia. • 1.47 

0.160 In) of ASTM-A7 (n-29) (2.06)* 

5. Topper, Sandor, and Morrow (13) ; 

strain controlled tests on 2024-T4 1.15 

aluminum (n«18) 

6 . Oowling (lU); 2024-T4 specimens; 
mostly In the elastic range, but 

some with large plastic strains; 0.363 

ra Inflow method used to count 
cycles (n*83) 


7. Topper and Sandor ( 15 ); strain 

controlled tests on 2024-T4 ^ 

aluminum with tensile and com- u-ujo 

pressive mean stresses (n«ll) 

3. Miner { 7); tests run at two or 

more stress levels on 245-T 0-080 

Alcad (n»ll) 

9. Schilling (16); welded cover plated 

beam specimens of A36 and ASTM AS 14 . ^ 

under simulated random bridge 
loading (n*36) 


1.30 0.962 0.627 

1.23 1.68 0.980 
1.23 0.341 0.269 
1-39 0.519 0.353 
1.14 0.186 0.161 

0.323 0.271 0.314 

0.309 0.219 0.262 
0.949 0.251 0.256 
1.26 0.78 0.54 


10. SAE fatigue Design A Evaluation 
Comnittee; tests on Man-Ten I 

RQC-lOO; notch specimen with 1.46 1.09 1.30 0.389 

cycl 1 c plasticity at notch root; 

See Table III; (n»54) 


Swanson flT^* amalgamation of random 
fatigue test data (n»671): Values 

based on an assumption that ■ 0.50, 1.21 1.04 3.726 3.50 

and F,(1.0) ■ 0.17 


♦Estimates of j, JSing a 'minimum SJI curve* 
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a keyhole notch. Failure was said to occur at crack initiation, which was 
arbitrarily defined as a crack of 2.54 inm (0.10 inch) A total of 54 
speciirans were tested. 

Results of several investigators in the 3AE program, using various 
analysis methods to predict cycles to failure, all using the PM rule, are 
presented in Table 6.2. Data from both materials are assumed to be from 
the same population. No attempt was made to judge the quality of the various 
analysis methods in computing the average values at the bottom of Table 6.2. 
These average values were used in Table 6.1. 

6.8.2 A Statistical Composite of A 

The goal of this study is to provide information to designers. To 
synthesize, for design purposes, the data of Table 6.1, statistics of com- 
posite data from two or more investigators, the following forms were used. 

, m 

A = S n.L (6.45) 

" i=l ' ' 

■ -ir ™ [(n, - 1)s,2 + n,(L - i)^] (6.46) 

A n-l 1 1 

where A is the sample mean of A, is the sample mean obtained by the ith 
investigator (as given in Table 6.1), n.j ii the sample size of the ith inves- 
tigator, m is the number of investigators, s^ is the sample variance of A, 
s^^ is the sample variance obtained by the ith investigator using A^ (this 
is what is given in Table 6.1), and n is the total number of observations by 
all investigators. Table 6.3 and Figure 6.9 provide summaries (using Equa- 
tions 6.45 and 6.46) of the analysis of the data listed in Table 6.1. 
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Tdble 6.2 

Sjrrr.ary Statistics an i Frcci 
SAE Cunulative Fatigue Damage 
Test Program (18) 


Socle and Mor-ow (igl 

Nominal strain life analysis 

Component calibration (incl 
mean) 

Neuber notch analysis 
Component calibration 

Dowling, Sroset Wilson ( 20 ) 

Base line case 

Base line case (Smith mean par.) 

Load strain extrapolation of 
elastic analysis 

Load strain calibration from 
Neuber 's rule 

Load, strain calibraiton from 
finite element analysis 

Constant amplitude load life 

Landqraf, Ricnards, La Pointe (21) 

Neuber notch analysis 
Load local strain analysis 
Strain based analysis 

Brose (22) 

Incl mean stress and overstrair 
effects 

AVERAGE VALUES 


t , 


--r 

A 






.385 

• 

.563 

1.06 

1.''9 

1.17 

.891 

1.00 

.855 

.90 

.693 

.745 

.823 

1.21 

.384 

1.14 

.939 

1.20 

.964 

.893 

.746 

1.15 

.922 

.848 

.739 

.687 

.505 

.632 

.920 

2.56 

2.C8 

1.32 

.713 

1.23 

.967 

.955 

. 779 

.677 

.561 

.457 

.674 

2.79 

2.22 

2.11 

.756 

2.69 

2.16 

2.01 

.748 

1.93 

1.51 

1.53 

.794 

1.43 

.970 

1 .56 

1.09 

1.46 

1 .09 

1 .30 

.389 
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A Statistical Suiimary of Models 
of Coiiiposite Data 
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DISTRIBUTION FUNCTION 
OF A,(%> 


^oincai Distribution Functions of Statistical 
Compo sites of Damage at Failure 
(Models are described in Table 6.31 


Composite 

AIuminum(4),y' 

/ //t . , 


/# 




I///'' 

i//0 Structural 
-'X/ Components (6? 

//'/ 

f /it 


DAMAGE AT FAILURE, A 



6.8 3 Analysis of the Variability in !\ 

For design purposes, it would be desirable to quantify the contribu- 
tions of various factors to the total variability in However, much 
of the detail regarding experimental procedure is lost in the data of 
Figure 6.6 and Table 6.1. Nevertheless, this data as a whole can be 
interpreted as an amalgamation of experiences and, as such, implicitly 
includes effects of many of the relevant factors. Values of coefficients 
of variation for material strength, e.g., yield strength of steel, ultimate 
strength of composites, etc. are typically less than 0.15 and, therefore, 
the values of listed in Tables 6.1 and 6.2 suggest a relatively large 
scatter in the data. 

Assuming the A is the product of several random factors whose dis- 
tributions are lognormal, it can be shown that 

1 + C.^ = n (1 + C.^) (6.47) 

^ i=l ^ 

where j is the number of factors which influence A, and is the coeffic- 
ient of variation of the ith factor. It is assumed that all factors are 
statistically independent. 

Often available is statistical information on fatigue variability inher- 
ent in the material, described by Cj^. Suppose that this is the only factor 
for which statistical information is available. Then Equation 6.47 reduces 
to 

1 + = (1 + (1 c/) (6.48) 

where C is the coefficient of variation of all of the other factors. 


6-25 


In general, it is difficult to separate the effect of mate’^ial 
variability from statistics on A using the method described above. In one 
example, the writer is suggesting that reasonable statistics for A when 
material variability is considered elsewhere is 

Median, A = 1 .00 

Cov, C, = 0.30 

with A being lognormal [8]. 

6.9 Statistical Considerations of Damage at 
Failure for Strainrange Partitioning 

Statistical analysis of damage at failure A for some strainrange parti- 
tioning (SRP) test data is presented in this section. This exercise was 
intended to provide necessary information for a total reliability analysis. 
Furthermore, statistics on A provide a quantitative measure of the perfor- 
mance of SRP. 

Fatigue damage after n cycles is (see Sect’’" 6.5 for definition of 
terms) 

4 

0 = n Z (f./N.) i = PP, CP, PC, CC (6.49) 

i=l ’ ^ 

Note that, in general, f^ and are random variables so that D also is a 
random variable. 

The SRP model defines a "predicted life" (cycles to failure) Np. In 

Equation 6.49 let D = 1 and n = Np at failure; choose the 

median, or best estimate, values f^ and . Then the predicted life is 

1 .Tf m ^ ■* = PP* ^P’ P^* 

TT " 

P 1 
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Consider a hypothetical experiment where several specimens are tested, 
having the same Np but different proportions of strainrange types f. , and 
therefore different lives . The cycles to failure for each specimen 
is recorded. The scatter in the results reflects the uncertainty in 
a) the life relationship, N • , i.e., scatter in Ae-N data, for each strain- 
range type, b) the process of computing the partitioned strainranges , f • , 
and c) in the performance of the linear damage model. Thus, is a 
random variable depending upon these three fundamental sources of uncer- 
tainty. 

Let damage at failure be a random variable denoted as This is 
the random variable which describes the uncertainty in the performance of 
the linear damage model. At failure D = A^ 3nd n = N^. An expression for 
can be written as (from Equation 6.49) 

^f ** V^(f-/N.) (6.51) 

11 1 

A f . and N. (i * U4) are random variables representing the three sources 
0 ’ i 1 

of uncertainty mentioned above. The distribution of A^ reflects only the 

inaccuracies of the SRP model. Unfortunately, there is no convenient way to 

measure A . However, it is possible from test data to obtain some informa- 
0 

tion on SRP performance. 

Equation 6.51 can be reformulated using the median (and constant) values 
f. and N.. 

Now A includes the variability of f- and N . , as well as model uncertainty. 
However, A can be measured directly from fatigue data. Substituting Equa- 
tion 6.50, it follows that 
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A = N^/Np 


(6.53) 


Given a single test specimen, the value of is 

*^i " ^0BS(i)/^'*PRE(i) 

where ^ = observed cycles to failure and = predicted cycles to 

failure. These values are available in the literature for SRP tests of various 
types. Table 6.4 summarizes the results of statistical analyses of several 
tests. For a given sample, A^., i = 1 ,r, the statistics were obtained using 
Equation 6.41 through 6.44. 

Following is a brief description of the differences in baseline and verifi- 
cation tests as noted in Table 6.4. 

Baseline Tests - Refer to the high temperature low cycle tests used directly in 
the establishment of the four SRP inelastic strainrange versus life - N) 

relationships. 

Verification Tests - Refers to the non-baseline high temperature low cycle fatigue 
tests and results used to check how well the established SRP relationships can be 
used to predict cyclic lives. A primary requisite of a verification test is that 
it should contain some feature or complexity not present in the baseline test. 

Finally it should be noted again that the scatter in A as reported by 
Table 6.4 includes uncertainties due to f^ and N. as well as uncertainties in the 
SRP model. It is generally desirable to separate these uncertainties ( to obtain 
the distribution of A^) in order to perform a reliability analysis. No attempt 
is made to do so in this study. 
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Table 6.4 


A Statistical Summary of SRP Data on Damage at Failure 


Sample 
Size, n 


Sample 

Mean 


Sample 

Std. 

Dev. 


Median 


AISI 316; 705C (1300F) 

Baseline tests 

Ref: Saltsman and Halford (23) 25 

AISI 316 

Verification tests 

Ref: Saltsman and Halford (23) 66 

AISI 304 

Verification tests 

Ref: Saltsman and Halford (23) 110 

AISI 304 and 316 

All verification test points 

Ref: Saltsman and Halford (23) 176 

AISI 316 

PP + PC + CC data only 

Ref: Saltsman and Halford (23) 9 

AISI 316 

PP + CP + CC data only 
Ref: Saltsman and 


AISI 316 

PP + CP data only 
Ref: Sa'tsman and 

AISI 316 

PP + PC data only 
Ref: Saltsman and 

AISI 304 

PP + CP data only 
Ref: Saltsman and 


Halford (23) 8 
Halford (23) 38 
Halford (23) 5 
Halford (23) 92 


1.25 1.11 .933 

2.87 5.81 1.27 

1.65 1.89 1.08 

3.02 6.96 1.21 

1.86 2.15 1 .21 

12.4 17.0 7.29 

2.67 5.09 1.24 

2.91 3.89 1.74 

1.71 2.03 1.10 


COV 

.891 

2.03 

1 .15 
2.30 

1 .16 
1 .37 
1.91 
1.34 
1.19 
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...Table 6.4 continued 


Mean 


Sample 
Size, n 


Std. 
Dev . 


Median COV 


AISI 304 

PP + PC data only 

Ref: Saltsman and Halford (23) 9 

lCr-lMo-.25V; 540C 

Normalized and tempered 

Ref’ Saltsman and Halford (23) 26 

lCr-lMO".25Vi 485C 

Normalized and tempered 

Ref: Saltsman and Halford (23) ^2 

2.25Cr-lMo; 540C 
Annealed 

Ref: Saltsman and Halford (23) 17 

2.25Cr-lMo; 540C 

Normalized and tempered 

Ref: Saltsman and Halford (23) 17 

2.25Cr-lMo; 485C 

Quenched and tempered 

Ref: Saltsman and Halford (23) 12 

AISI 304; 650C 
Solution annealed 

Ref: Saltsman and Halford (23) 15 

AISI 304; 565C 
Solution annealed 

Ref: Saltsman and Halford (23) 13 

Rene 95; 1200F 
Baseline Tests 

Ref: Hyzak and Bernstein (24) ^2 

Rene 95 

Verification Tests 

Ref: Hyzak and Bernstein (24) 13 


1.54 1.16 

1.00 .409 

1.27 .293 

1.21 .241 

0.820 .190 

0.615 .222 

0.899 .367 

0.467 .185 

1.25 .830 

0.920 .479 


1.23 .750 

.928 .408 

1 .24 .231 

1.19 .198 

.799 .232 

.579 .350 

.832 .408 

.436 .395 

1.04 .663 

.816 .520 


16 


Rene 95 

Validation Tests 

Ref: Hyzak and Bernstein ( 24 ) 
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1.13 


.327 1.08 


.290 


6.10 Statistical Distribution of ^ 


The least squares program for distribution analysis as described in Section 
2.5 was used to analyze the data on A. The results are summarized in Table 6.5. 

The model which best fits the data for each sample is indicated by the box 
in Table 6.5. The results show clearly that the lognormal is generally the best 
fit model for the samples considered. It should be noted that this implies that 
it is reasonable to assume that c-ycles to failure is lognomially distributed. 

6.11 Summary Comments on Chapter 6 

Provided in this chapter were equations, using the concept of damage, 
which are used by designers to predict fatigue failure under both variable 
and constant amplitude stresses. Statistical information on the distribution 
of damage at failure was also presented. Summary statistics provide a quan- 
titative description of the performance of Miner's rule. Statistics were 
also provided to describe the behavior of the SRP model. 

Given that uncertainty exists in all of the fatigue design factors, 
the question now is how to use the damage equations to make decisions which 
will produce a high reliability product. Chapter 7 will sumnarize elements 
of modern structural and mechanical reliability theory, and Chapter 8 will 
address specifically the problem cf fatigue reliability. 
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Table 6.5 


A Summary of the Best Fit Distrioution for SRP Data on Damage 
at Failure Using the Lea it Squares Method 


Samole 
Size, n 


NORMAL 


AISI 316; 705C (1300F) 

Baseline tests 

Ref: Saltsman and Halford (23) 25 

AISI 316 

Verification tests 

Ref: Saltsman and Halford (23) 66 

AISI 304 

Verification tests 

Ref: Saltsman and Halford (23) 110 

AISI 304 and 316 

All verification test points 

Ref: Saltsman and Halford (23) 176 

AISI 316 

PP + PC + CC data only 

Ref: Saltsman and Halford (23) 9 

AISI 316 

PP + CP CC data only 

Ref: Saltsman and Halford (23) 3 

AISI 316 ^ 

PP + CP data only 

Ref: Saltsman and Halford (23) 38 

AISI 316 

PP + PC data only 

Ref: Saltsman and Halford (23) 5 

AISI 304 

PP + CP data only 

Ref; Saltsman and Halford (23) 92 


.551 

.731 

.553 

.772 

.484 

.485 

.715 

.578 

.566 


- A 

e * /I - D 

LOGNORMAL 
j .197| 


.108 


.124 


.147 


WEIBULL 

.280 

.222 

.178 

.237 


.179 
.332 
.141 I 


.282 



I .146 
.320 
.251 
.370 
.204 


* smallest value of c indicates that distribution which best fits the data; 
the box indicates the best fit 
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...Table 6.5 continued 


Sample formal LOGNORMAL 
Size, n 


AISI 30*1 

PP + PC data only 

Ref: Saltsman and Halford (23) 9 

lCr-lMo-.25V; 540C 

Normalized and tempered 

Ref: Saltsman and Halford (23) 26 

lCr-lMo-.25V; 485C 

Normalized and tempered 

Ref: Saltsman and Halford (23) 12 

2.25Cr-lMo; 540C 
Annealed 

Ref: Saltsman and Halford (23) 17 

2.25Cr-lMo; 540C 

Normalized and tempered 

Ref: Saltsman and Halford (23) 17 

2.25Cr-lMo; 485C 

Quenched and tempered 

Ref: Saltsman and Halford (23) 12 

AISI 304; 650C 
Solution annealed 

Ref: Saltsman and Halford (23) 15 

AISI 304; 565C 
Solution annealed 

Ref: Saltsman and Halford (23) 13 

Rene 95; 1200F 
Baseline Tests 

Ref: Hyzak and Bernstein (24) 42 

Rene 95 

Verification Tests 

Ref: Hyzak and Bernstein (24) 13 

Rene 95 

Val idation Tests 

Ref: Hyzak and Bernstein (24) 16 


I .304 
.313 


.423 


.228 
I 


.260 


.247 


.167 


.417 

.285 

.245 


.272 

.420 

.221 

.622 


.231 


i .265 

» 


I .212 


i .286 


.218 


.414 


.171 , .109 I 


WEIBULL 

.365 

.257 

.296 

.291 

.188 

.238 

.434 

.269 

.421 

.268 

.191 
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Chapter 7 A REVIEW OF MODERN PROBABILISTIC DESIGN THEORY 


7.1 Introduction 

An example cf a fatigue reliability problem is provided in Table 7.1. The 
design factors are treated as random variables reflecting uncertainties due to 
statistical scatter in observations as well as assumptions made in the analysis. 
The probability of failure is computed by solving a probability problem. If p^ 
is less than p^, a predetermined maximum allowable risk, the design is considered 
to be safe. 

It is assumed that the condition for failure is described by an algebraic 
expression. For example, from Table 7.1, the condition for fatigue failure is 

YyKSQL/bh^)]"'" < N^ 

Each term on the left hand side could be treated as a random variable. The gen- 
eral goal of design is to select design variables so that the probability of 
occurrence e.g. the above event, is acceptably small over the design life. 

It is the goal of this chapter to provide a concise summary of the tools 
available for reliability analysis, i.e. computing the reliability or the proba- 
bility of failure of an index thereof. A number of methods are reviewed. No 
specific recommendation is made regarding a preferred approach. Indeed a relia- 
bility analysis must be tailored to the design procedure, the data available, and 
the way in which the results are to be used. 

The following summary of modern probabilistic design theory is based on the 
articles by Hasofer and Lind (1), Ditlevsen and Skov (2), and the reports of 
El 1 ingwood, • Galambos, MacGregor, and Cornell (3), Comite Europeen du Beton No. 112 
(particularly the commentary by Rackwitz) (4), and CIRIA No. 63 (5). Many of the 
operations described are application of elementary probability theory. Therefore 
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Table 7.1 


A SUMMARY OF A FATIGUE RELIABILITY PROBLEM 



Stress Amplitude, ^ ® 


(b) Fatigue Strength of Material ;NS"’ = K 
where N = number of cycles to failure 

(c) Design Life, 

(d) Design factors (Q, L, b, h, K, m) 
are in general, random variables. 

(e) Non-statistical uncertainty exists 

i Due to the effects of manufacturing and assembly operations and use in 
service 
K = YK 

0 

ii Due to assumptions made to compute fatigue stresses 
S = ZS 

° Y and Z are random variables which describe this 'professional" 
uncertainty 

FIND The probability of failure, p^ 

formulation The event of failure is N < The probability of failure is 

Pf = P[N < N^] 

Substituting into the above expression for N, 

Pf = P{YK [Z(6QL/bh2)]“"’ < N } 

We are left with what is, in general, a difficult probability probleir 
Moreover we must also define a safe level of risk p . The design is 
acceptable if Pr < p ° 

‘ ■ o‘ 


7-2 



the reader is encouraged to have available a good elementary text in probability 
theory and statistics, e.g. Meyer (6), Hines and Montgomery (7), Bowker and Lie- 
berman (8), etc. 

7.2 Basic Definitions 

A basic assumption in the following discussions is that the vailure condition 
can be written as an algebraic function, 

R(X) < S(Y) (7.1) 

where R and S, functions of several random design factors X and Y, represent 
strength and stress functions respectively. The assignment of factors to R and 
S is often arbitrary. The tilde underneath the variable indicates that the quan- 
tity is a vector (or multivalued). The design factors U * (^*7) can be separated 
into the known random variables V and the design parameters to be determined A, 
e.g., cross sectional area, sectional modulus, etc. Thus U = (V»^)* Following 
are the basic assumptions regarding U. 

It will be assumed for simplicity that the standard deviation = 0 for each 
A^. Geometric variances are usually small compared to others. Furthermore, it is 
assumed that all U. are independent. If two or more design factors are correlated, 
then a transformation to uncorrelated variables may be made (1). 

Define the fai lure surface or 1 imi t state as 

g(U) = R(X) - S(Y) = 0 (7.2) 

The failure surface is the boundary between the safe and the failed regions in 
design parameter space. In general, there can be several limit states for a com- 
ponent, e.g., yield, buckling, fatigue, etc. 
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The failure function , g, is defined as 


g(y) = R(X) - S(p (7,3) 

condition (or failure event ) occurs when S(Y) exceeds R(X). Thus 
the failure condition is, 

f 0 (7.3) 

For example. Figure 7.1 shows a single element, single load structure with 
random design factors R and S. Here R and S are considered the basic variables. 
The overlap in pdf's suggests that it is possible to select R and S at random 
such that R < S and failure occurs. 

Design parmeter space for this case is shown in Figure 7.2. The limit state 
is R » S. The failure condition is R - S < 0, or simply R < S. 

THe Full Distributional Approach 

It IS assumed that the exact statistical distribution of each design factor 
U^. is known. 

mipbability of failure (or ri^) p^ is defined as the probability of the 
failure event, 

Pf = P(g(U) < 0 ) (7,4) 

For example, ,f R and S are the basic variables, the probability of failure is 
P^ = P(R - S < 0) or P(R < S). 

In general, p^ can be computed as (3) 

^f ••• ••• ^*^k (7.5) 

where is the region of U where g(U) < 0, f^j is the joint pdf of U, and k is 


A, ,?rea 

R, strength 

S, stress 
(Q/A) 



Figure 7.1 Probability Density Function? (pdf) of 
Stress and Strength 



Figure 7.2 Design Parar.eter Space when R and S are the Basic Variables. 




the number of design factors. In the special case where the two basic variables 
are (R,S) and statistically independent, it can be shown that 

00 

Pf = 1 - r f|^(r)Fj(r) dr (7-6) 

= f fj(s)Fj^(s) ds 

0 

where f(-) = probability density function and F(*) = cumulative distribution 
function. The subscripts denote the random variables. 

Examples CZo6ed FoAm SoZu.tiom 

(a) Only one random variable . If S = S^, a deterministic real number, then 
Equation 7.6 becomes, 

Pf = P(R < Sq) = F^(Sq) (7.7) 

and if R = R^, deterministic with only S random, 

Pf = P(Rq < S) = 1 - FgCR^) (7.8) 

Example: Cycles to failure N of a component is known to have a Weibull dis- 
tribution with parameters a = 3.58 and 6 = 21.25 thou-cycles. N is a "strength" 
variable. The design life ("stress") is given as = 12 thou-cycles. Compute 
the probability of failure, i.e. if a member is selected at random, what is the 
probability that its life will be less than 12 thou-cycles? 

Pf = (N < Nq) = Fp(Nq) (7.9) 

= 1 - exp[12/21.25)^-^®] 

= 0.121 

The pdf of N and p is illustrated in the following figure. 
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( 5 ) The Normal Format ; Assume all U^. have normal distributions and that the 
limit state has the linear form. Let 


k 

Z = g(U) = A +.1^3. U. 


(7.10) 


where A and are constants. 

From elementary probability theory (e.g. Reference 7) the random vari.ible Z 
will also have a normal distribution. The mean and standard deviation of Z are 


k 

y, = A + [ B.u. 
^ i=l ^ ^ 


a 


Z 




2 


(7.11) 


The failure condition is Z < 0. The probability of failure is 


= P(Z < 0) 


(7.12) 


The probability density function of Z is shown in Figure 7,3. Using stan- 
dard methods of calculating probabilities for a normal variate. 


P^ = P(Z < 0) 


a - 


u 




z < n 

a 


ZJ 


(7.13) 
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where <I> is the standard normal cdf, and 

6 = (7.14) 

3 is called the s afety index . 6 is related to as shown in Figure 7.4. Note 
from both Figure 7.3 and 7.4 that the safety index 6 plays the same role in pro- 
babilistic design as does the factor of safety in conventional design, i.e., the 
larger is 8, the smaller is the risk, or p^. 

Example of the Normal Format: Consider the tensile bar of Figure 7.1. 

Stress, S * Q/.^ where Q is a random load variable having mean Uq and standard 
deviation Oq. The strength of the bar R is considered as a random variable having 
mean and standard deviation of and Op. Compute the probability of failure. 

The mean and standard deviation of Z is, (from Equation 7.11) 

Uz . - Pg/A 

Then, 

= $(-B) (7.17) 


(7.15) 

(7.16) 


where 8 = 

(c) The Lognormal Fonnat . The lognormal distribution plays an important role 
in probabilistic design and in particular for problems which involve fatigue. 
Properties of the lognormal distribution are summarized in Table 2.5 with a more 
complete discussion in Appendix 1. 

Consider a restatement of the limit state and a redefinition of g. Let, 
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g(u) = R(X)/S(Y) = 1 


(7.18) 


Now the failure condition is g < 1. Assume that (a) all have lognormal dis- 
tributions, and (b) g(U) is a multiplicative function of the design factors 


g(y) U.‘ 


(7.19) 


where B and all a. are constants. Let Z = In g. 


Z 


K 

In B + y 
i=l 


a^ln 


(7.20) 


Because U. is lognormal, it follows that InU^is normal, 
comes identical to the normal format. 

The failure condition is g < 1 or, as above, Z < 0- 
Equation 7.17 with B * u-r/a^, and 


u = Z = In g = U.^‘ 


Then this format be 


The is given by 


(7.21) 


In 


n (1 + C.^)^i 
i=l 1 


where the tilde indicates a median value and C^. is the coefficient of variation 
(cov) of U^. . For a random variable X, the median and cov in terms of the mean 
and standard deviation are respectively. 




X = uy/Ji + c 


= 


(7.22) 


See Appendix 1 for additional detail. 

The lognormal format is particularly useful in fatigue design because some 
of the important equations tend to be multiplicative. 
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Example: The limit state for fatigue under a stationary gaussian process having 
RMS, a is (see Equation 6.29) 

A = ^ (7.23) 

where A = N^(/?)^r(m/2 + 1), a constant. A, K and a are random variables. Fail- 
ure occurs when the right hand side (stress) exceeds the left hand side (strength) 

g(U)=A!L (7.24) 

Aa 

failure occurs when g < !• The mean value of Z is, 

* ln(AK/Aa"’) (7.25) 

where the tilde indicates median values. Note that for th.' lognorm al, the m ean 

Uyand the median x of the random variable X are related by X = 

Also » 2 , 

0 ^ = [ln{(l + C^^)(l + C^^)(l + (7.26) 

Then 3 and p^ are determined by Equation 7.13 and 7.14. 

7.4 The Full Distributional Approach: Monte Carlo Methods 

Consider a complicated design problem having a limit state, for example, 

Z = g(U) = (AX®/Y^) + CU^lnV + K = 0 (7.27) 

in which A, C, a, b, c, and K are constants and X, Y, U and V are random variables 
having a Weibull, normal, lognormal and gamma distribution respectively. Evalua- 
tion of p^ using Equation 7.5 is a practical impossibility. 

Monte Carlo methods, however, can be used to obtain approximate solutions to 
complicated probabilistic design problems as described above. Estimates of the 


'* 

=i 
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distributions of complicated functions of random variables can be used to evaluate 
risk. 

Monte Carlo methods are widely used, and only a demonstration will be pre- 
sented as follows: 

Example: Cycles to failure in low cycle fatigue for AISI 316 stainless steel 
is given as 

N = (7.28) 

♦ 

where Zi, G and c are lognormal ly distributed random variables having the following 
values for the median and coefficient of variation. 



Median 

COV 

A 

1.00 

0.30 

G 

0.222 

0.593 

£ 

0.015 

0.30 


The design life is given as = 40 cycles. Determine the probability of failure. 

The limit state is N = N^, and the event of failure is (N < N^) . Thus 

Pf = P(N < (7.29) 

An estimate of the cdf of N can be made by first sampling K values of A, G 
and e at random from the lognormal. These values are used to compute a sample 
of size J of N from Equation 7.28. Statistical analysis of N can include esta- 
blishing the empirical cdf. For the above example the empirical cdf is automa- 
tically plotted as shown in Figure 7.5. The plot shows an estimate of p^ = 0.0054. 

A closed form solution to the above problem is available (See Appendix 1 ). 
However, if for example A had a lognormal, K had a Weibull and e had a normal, 
establishing the distribution of N by analytical methods would be a practical 
impossibility. Nevertheless, the approximate Monte Carlo solution is no more 
difficult than that of the above example. 
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Figure 7.5 

The Empirical Distribution of N Plotted on Lognormal Probability Paper 
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7.5 First Order Second Moment Methods 


In the full distributional approach described above, it is assumed that each 
design factor has a known distribution. However in practice, these probability 
laws are seldom known precisely because of a general scarcity of data. In some 
cases only the mean and standard deviation may be known with any confidence and 
even if the distributions of all design factors are known, it is usually imprac- 
tical to evaluate Equation 7.5. 

The difficulties above have motivated the development of first order second 
moment (FOSM) reliability methods. The random variables are characterized by 
their mean and standard deviation. 

The basic measure of reliability is the safety index B = exam- 

ination of Figure 7.6 it is clear that B as a measure of safety can be employed 
for any g(U) and any distribution of each U^. . However Z has a normal distribution 
and Equation 7.13 for p^ applies only for the linear form. Equation 7.10. 

Determination of and where Z = g(U), a complicated function, is in 
general very difficult. However these terms are relatively easy to evaluate if 
g(U) can be linearized by a Taylor's series expansion about U . Excluding higher 
order terms , n 

I = g(U*) + I{'ag/3U.),*(U. - U *) (7.30) 

i=l M ^ ^ 


It is now a relatively simple matter to compute 3nd a^, but a key consideration 
is selection of the appropriate linearizing point U . 

Mean l/o£ue MeC/iodi. In earlier structural reliability studies, the point U* 
was set equal to the mean value y. Assuming that all U. are independent, it can 
be shown that 


a 


Z 



2 

U 1 


(7.31) 
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Figure 7.6 


Relationship with Safety Index to Risk and to u, and j, 

L L 

Failure Z = g(U) < 0 
Mean and Standard Deviation of Z; 

'>^1 ~ E[9(U)] ^ 2 , ~ /var[g(U) ] 


Safety Index; 

3 




Changes in (larger) and 
(smaller) to increase S, will 
decrease p^r. ••• for any distri 
bution of Z. 



where = the standard deviation of U^. . 

Example: Consider the tension bar of Figure 7.7. Compute the safety index. 

Z - g(U) = g(R,A) = R - Q/A (7.32) 

Note that the function g is formulated so that the limit state is Z = 0 and the 
failure condition is Z < 0. 




(7.33) 


rr 2 - 2 . 2 


with the partial s evaluated at the means. 




(7.34) 


Substituting the values from Figure 7.7 u^ = 26.3 and Oj ~ 6.45. Thus 


S = = 4.07 (7.35) 

Evaluation of 6 for this case was relatively easy. Even for more complicated 
limit state functions the operations are not Jifficult. On the other hand there 
are significant limitations and criticism of mean value first order second moment 
methods which have focused on the following: 

(a) In this approach, the measure of safety is 6, which in some cases gives 
weak information on the probability of failure, 

(b) Information on distributions, if available, cannot be included in a 
logical way, 

(c) The linear approximation of g(U) at the mean of the various variables 
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Figure 7.7 

data for example problem 


STRENGTH 5 R 
* 62 ksi 

» 6.2 ksi 
AREA, A 

. 2 

* 2.80 in 

* 0.14 in^ (very poor quality control) 


1 

LOAD, Q = 100 kips (Deterministic) 

Limit State: 

Formulation 1 

Z = R - Q/A * 0 
Formulation 2 

Z « RA - Q * 0 
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appears to be too inaccurate in the face of severe non-linearities of this 
function in many design situations, i.e., the higher order terms in the Tay- 
lor's series expansion are important. 

(d) Equation 7.31 for ^nd a, is valid only for relatively small variances 
(typically for C^. < 0.15). In the fatigue problem most design factors have 
"large" variances, ... again the higher order terms are important 

(e) and the most important, the method fails to be invariant in that S depends 
upon the mechanical formulation of the problem. 

The lack of invariance problem can be illustrated using the above example. 
Note that an equally valid statement of the failure condition is 

Z = RA - Q < 0 (7.36) 

Using procedures described above. 


- Q 


(7.37) 


2 2 2 ^ 22 
h - Ua 


Substituting values from Figure 7.7, 73.6 and = 19.4 and 6 = 3,79. 

Different values of 6 are obtained from two equally valid statements of the 

failure condition. It is principally this problem with the mean value method 

that led to the development of the generalized safety index (1). 

Advanced f<A6t OideA Second Moment Methods. It is now recognized that S will 

* 

be invariant to formulation if the linearization point U is chosen as a point on 

the failure surface. The following discussion describes how to compute the genera- 

★ 

lized safety index and to establish the appropriate U . 

Hasofer and Lind have proposed a general approach which is exact and invariant 
to mechanical formulation (1). As an example, consider R and S as the basic 
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variables . In Figure 7.8 the probability of failure will be the volume under 
the joint density function in region n. Perform a transformation ; ds^lne the 
reduced variables : 

r = (R - S = (S - U 5)/'^5 (7.38) 

r and s have means equal co zero and standard deviations equal to one. 

Figure 7.9 shows the space of reduced variables. Note that again p^ is the 

volume under the joint density function in n. Thus the distance A = OP is a 

measure of reliability. As the failure surface moves to the right, p^ gets smaller 

and A gets larger. Therefore A is a measure of reliability. Furthermore from 

geometric considerations 

A ■ (ur - 

which is exactly the same as 8 obtained from Equation 7.31 for g = R - S. 

The design point is defined as the point on the failure surface closest to 
the origin of the reduced coordinates. In the above example it would be the 
point (s*,r*) as shown in Figure 7.9. The design point in the basic variables is 

* * X 

then determined by inverting the original expression, e.g., R = r from 

Equation 7.38. 

Define the generalized safety index , 3 as the minimum distance from the 
origin of the reduced variables to the failure surface 

/ 2 2 2 

U] + U2 + ...u^ } (7.40) 

Subject to g(U) = 0, where the reduced variables in terms of are 

u^ = (U. - U^)/a. (7.41) 
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Fig 


Parameter Space When R and S are the Basic Variables and the Limit 

State is Z = R - S = 0 



R 


Contours of joint pdf of R & S, f 







MICROCOPY RESOLUTION TEST CHART 
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Such a procedure will produce a design algorithm which is invariant to formulation 
because the failure surface is the same regardless of hew it was written. Figure 
7.10 illustrates the failure surface in two dimensions. U. are the '•educed varia- 
bles. 

Another viewpoint of this problem is to consider the system of equations (6), 


3g^/3u^ 


(7.42) 


u, = -a,B 


g-|(u.j ) = 0 

where 8 is the distance to point u on the failure surface g-j in reduced woordi- 
nates. is the direction cosine of the i‘" coordinate relative to u*. These 
systems of equations are so'ved for the o/s and u* which minimizes 3. 

It is also possible to redefine any or all of the basic variables by making 
a log transformation U! = In U.. Then, the mean and standard deviation of U' 
are (assuming the lognormal distribution moments) 


Uy,- In uy - 1/2 In (1 r Cy^) = In U 


(7.43) 


oy, = /in ( 1 f Cy2) 


For convenience, the subscript "i“ has been dropped. The reduced variables are 


u = (u'- Viy»)/c^y« 


(7.44) 


so that the original variables in terms of the reduced variables are 
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Figure 7.10 

Generalized Safety Index in Two Dimensions 





U = exp(;;yi+ uJyt) 


(7.45) 


In general the log transformation would be performed on large variance 
design factors, and in particular would be useful in the fatigue problem. 

To compute one must in general perform an optimization calculation. The 
author's experience is that most failure surfaces are well behaved and as a 
result, convergence to an accurate (four place) result using a relatively crude 
optimization algorithm, reguires a negligable amount of computer time, even for 
several design variables. 

Example: Consider the problem of Figure 7.7, The two valid expressions for 
the failure surface are 

Z = R - Q/A = 0 (7.46) 

» RA - Q * 0 


The reduced coordinates are 

r = (R - UR)/Of^ a = (A - u^)/a^ (7.47) 

Using the values in Figure 7.7 and substituting for R and A into either of the 
expressions for Z, the same expression for the limit state results 

g.|(r,s) = (62 + 6.2r)(2.8 + .14a) - 100 = 0 (7.48) 

This function is plotted in Figure 7.11. Also shown is 8, the minimum distance 
to the failure surface as well as the design point (a , r ). Here the values were 
determined graphically. 

Note that from Figure 7.10, the procedure for determination of 3 is equiva ■ 

★ 

lent to linearizing the limit state equation at u and computing 3 associated with 
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the linearized rather than the original limit state. From Equation 7.30, 
it can be shown that 

=2^ = E[(Z - = I OZ/3U,)'*0j^ (7.50) 

★ ★ 

In the above example of Figure 7.11 R = 38.0, A = 2.632 and 8 = 4.05 in 
both formulations. 

In summary, the generalized safety index 8 provides a measure of reliability 
for which is invariant to the mechanical formation of the limit state. It is 
necessary to specify only the mean and standard deviation of each design variable. 
The S so defined contains in general, weak probability information. However, if 
the variables are normal or lognormal (and the log transformation is used), then 
~ <^(-6) is a good approximation to the exact probability of failure as would 
be computed by Equation 7.5. 

7.6 Extended Form of the Generalized Safety Index 

The procedure described in the previous section gives values of 8 which may 
be related to the probability of failure only in those special cases where design 
factors have normal or lognormal distributions. But many design problems involve 
random variables having other distributions e.g. Weibull. Computation of the 
probability of failure in the general case of a non-linear limit state and design 
variables having arbitrary distributions can be accomplished by a very complicated 
numerical integration (Equation 7.5) or by Monte Carlo methods. Neither method 
is practical for general design application. 

In 1976 Rackwitz has proposed an extension of the generalized safety index 
concept (6,7,7,12) in which appropriate distributional information could be 
incorporated. The basic idea is to transform the non-normal variables into 


7-27 


equivalent normal variables. This transformation may be accomplished by approx 


imating the true distribution of variable U^. by a normal distribution at the 
★ 

value corresponding to a point on the failure surface. The justification 
for this is that if the normalization takes place at the point close to that 
where failure is most likely, (i.e., minimum 3 ), the estimates of the failure 
probability obtained by the approximate procedure should approximate the true 
(but unknown) failure probability quite closely (6). 

The mean and standard deviation of the equivalent normal variable are 
determined such that at the value U^t the cumulative probability and probability 
density of the actual and approximating normal variable are equal. Thus, 

, ct(| ^(F.(U.)]) (2.24a) 

u'J . U* • f ^[F.(U*)]aJ (2.24b) 


in which F^ and f^ * non-normal distribution and density function of and 
p(-) is the density function for the standard normal variate. Having deter- 

N N 

mined and o^ of the equivalent normal distributions, the solution proceeds 

as described in Equations 7.40 to 7.42. Inasmuch as the checking point 

* N N 

variable changes with each iteration, the parameters u^. and must be 

recomputed during each iteration cycle also. However, since all calculations 
are performed by computer, this does not materially add to the complexity of 
the reliability analysis described earlier. 

It has been observed that this approximate technique often yields excel- 
lent agreement with the exact solution, e.g. as obtained by Equation 7.5 (6.12). 
However, it has been noted (15) that the checking point may not correspond 
exactly to the point where the joint probability density is maximum and failure 
is most likely. However, this procedure does not reduce the error which is due 
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to the linearization of what may be a generally nonlinear failure boundary 
at the design point. Unless the failure boundary is highly nonlinear, how- 
ever, as is the case in some stability problems, this source of error is 
small compared to the accuracy with which most of the parameters in engineer- 
ing reliability analysis can be estimated. 

A computer code using the Rackwitz algorithm has been written as an 
activity of this study (called SUPER-HASL). An example which illustrates 
the efficiency of this algorithm is presented in Section 8.4. In this demon- 
stration, the probability of failure is computed at roughly 1/6 of the cost 
of Monte Carlo. 

SUPER-HASL has only recently been developed. The efficiency and accuracy 
of the Rackwitz algorithm will be investigated in future studies. 

7.7 Summary Comments on Target Reliability Levels and Technique 

Important questions have not been addressed in the previous discussion. Fol- 
lowing are answers to some key questions regarding implementation of reliability 
methods . 

Hotc ,V!ach Con^id&ncz Can We Put htta p^? 

A computed value of probability of failure p^ should be interpreted as an 
"index" of risk rather than an actuarial prediction, p^ can be sensitive to the 
many assumptions used that one may have little confidence in the quality of the 
point estimate of p^. An estimate of p^ is sometimes called the "notional pro- 
bability of failure". In practice, 8 is now more commonly used (See Ref 4). 

What Should be a Tatget Value, 6 (o-t p^)? 

The target design value of 3, denoted as 3^, depends entire y on the appli- 
cation. Values of 3 implied in current codes and procedures typ cally range from 
2.0 to 6.0. The value of 8 = 3.0 has been suggested for general construction (4). 
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To establish tne engineer can (a) attempt to evaluate S inherent in current 
design procedures as a reference, and (b) use an implied relationship betv/een B 
and p^ assuming that p^ accurately describes risk. Ultimately the designer has to 
use his judgement and experience in making a decision regarding acceptable levels. 

Can Wc Provide a Gzncxal Pztiab-^ ^itu Met/icd Th-’t Can be Used Fo'i AiZ Cases? 

It is impossible to recommend a general procedure which can be used in all 

cases because each class of problem seems to have its own unique characteristics. 
A reliability method must be tailored to the sit ation. However some general 
observations can be made from the above discussion. A summary of procedures 
available for each problem type is provided in Table 7.2. 
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Table 7.2 


A SUMMARY OF METHODS OF PROBABILISTIC DESIGN 


Form of the Failure Function 
Z = g(u) 

Reliability Method of Preference 

A stress and strength variable have 
known distributions 
Linear (Equation 7.10); all design 
factors normal 

can be computed by numerical integratiO' 
^(Equation 7.5) 

Z is exactly normal and p^ can be computed 
exactly (Equation 7.13 anU 7.14) 

Multiplicative (Equation 7.19); 
all design factors lognormal 

Z ® In g is exactly normal and p. can be 
computed exactly (Equation 7.13 and 7.14) 

1 

g is any function; all design factors 
have known distributions 

Monte Carlo method; a relatively easy 
method for estimating p, but unweildly for 
a designer; the Rackwitz algorithm may be tli 

g is any function; all design factors 
have known means and variances but 
the variance is small (COV < 0.15) 

besl 

The mean value first order se-ond moment 
method provides an easy approximate 
method to evaluate S. 

g is any function; all design factors 
have known means and variances; if 
variances are large, use a log trans- 
formation on each factor 

The generalized safety index can be 
computed using a digital computer code 
(Ref 2,6). 

g is any function; some design factors 
have known distributions but others 
have only means and variances known. 

The generalized safety index can be com- 
puted. A digital computer program is 
required. The method of analysis is not 
described herein, but is provided in 
Ref. 6. 
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Chapter 8 DEMONSTRATION OF RELIABILITY ANALYSIS TO FATIGUE PROBLEMS 
8.1 Preliminary Remarks 

Each design and safety check problem seems to have its own unique features. 
Therefore it would be difficult to propose a general reliability format which 
could be used in all problems. 

Three examples, presented in this chapter, illustrate how some of the relia- 
bility methods discussed in Chapter 7 can be employed in fatigue problems. The 
examples treat models used to predict high temperature low cycle fatigue. The 
first uses strainrange partitioning to define the limit state and the second 

and third use the local strain approach. 

These examples do not exhaust available reliability methods but provide only 

a sample of the possibilities. 


8.2 Example 1; Strainrange Partitioning in a Reliability Format . 

Consider a component for which high temperature fatigue analysis is performed 
using strainrange partitioning. Only PP and PC strains are present. The limit 
state is 


A = N, 


PP 


+ ^PC 




( 8 . 1 ) 


See Section 6.5 and 6.9 and Reference 1 for definition of terms. The data for 
all parameters are given in Table 8.1. 

The right hand side of Equation 8.1 is damage at the design life N^. Fail- 
ure occurs if this exceeds A, the damage at failure. Thus the left hand side can 
be considered as "strength" and the right hand side as "stress". 

The random variables A, G, H, and Ac are related in such a way that a closed 
form solution for probability of failure, e.g. using Equation 7.5, would be 
extremely difficult. In lieu of this, two approximate techniques for reliability 


analysis are pursued. 
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Table 8.1 

Data for Example 1 


500 cycles, design life 
0.70 
0.30 
-1 .711 
-1 .188 

Random Variables (all assumed to be lognormally distributed) 


Median | 

Coefficient of Variation 

A 

1.0 

0.30 

G 

.222 

0.406 

H 

1.673 

0.393 

Ae 

0.02 

0.300 


PP 


PC 


Y 

n 
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1. The generalized safety index can be computed. A log transformation is per- 

★ 

formed on all four variables. A digital computer code (HASL) is employed. 

The output is presented in Table 8.2. The computed safety index is b = 4.08. 

The generalized safety index provides a good approximation to the probability 
of failure using = '5(-6) when the design factors are normal or log- 
normal, even though the limit state is not linear in the design variables 

Here p^ = 2.28 x 10'^. (See Figure 7.4) 

2. The Monte Carlo method can be used to estimate the probability of failure. 

Equation 8.1 can be rewritten as 


N = A/ 


PP 


PC 


3(Ac)^ H(Ae)"J 


= N. 


( 8 . 2 ) 


N is a random variable denoting cycles to failure. Failure occurs if N < N^. 
The Monte Carlo method (Section 7.3) is used to obtain a random sample of N, 
the distribution function of which is plotted on lognormal probability paper 
in Figure 8.1 . 

Graphically extrapolating the empirical distribution function to = 500 
cycles, the inverse cdf, equal to -3, indicates that 3 == 3.9. The corresponding 
p^ = $(-3.9) = 4.85 x 10'^. 

Considering the degree of extrapolation required for the Monte Carlo method, 
it is not surprising that there is only fair agreement between the two methods. 


*See Section 7.4 for the description of the generalized safety index. 

**Just how good this approximation is, in general, is presently under study. 
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Table 8.2 


Output of HASL 
Safety 


A Program for Computing the Generalized 
Index--for the Example of Table 8.1. 


DESIGN VARIABLES 


VARIABLE 

transformation^*) 

NEAN/H60IAN 

STo oev/cov (<^) 

...DELTA 

i 

. .I0p00f4-01 

.30QOOE>00.. 

G 

I 

• 22200E«>00 

•A0600E400 

— S.PlkN 

I 

— ti5QQj3JE-02 


H 

1 

.16730E^01 

— f J V. W \d \J c ▼ w w 

' .39300E^00 


DESIGN 

POINT 


V AK 1 

ffEOUTCny VALUE 

TaSTC-VaTOIE 

DELTA 

XRm* -2.02318 

X( !)• 

.95215E^00 


KITTD- -l,T92n — 

xm- 

.i"2J9V£ + 0O 

EOSLN 

' M. 

XR(3»« 2.99219 

X ' 3)- 

..48140E-02 


XRT^l* -IT16406 

X(4T.' 

■ ;T(Tr$Te‘>oi 


SAFETY INDEX# beta • 4,07832 


NOTES: 

(a) HASL is a digital computer program developed by P. H. Wirsching, University 
of Arizona, for computing the generalized safety index 

(b) A "1" indicates a log transformation on each variable 

(c) /li th a log transformation, the listed values are the Median and the COV 
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INVERSE cor or N C-p; 


.999 
.998 

- . 995 

-.99 
L.93 


.95 

..90 

t .30 
.70 
.60 
. .50 
j-.40 
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o 

\ 




o 

« 


.20 

..10 

.05 


t .02 

.01 
.005 

. . 002 ' 

- .001 
L .0005 


/ 

/ 

/ _ 

’ \ 


Figure 8.1 

Enpirical Distribution Function of Cycles to Failure 
N Plotted on Lognormal Probability Paper 


-3.9 


10 


y 


/ 






/ 


y 






/ 


‘ = $(-3.9) 

= 4.85 X 10' 


-I r 


CYCLES F^.^'J^E, *! 


in* 


n = 500 


Extrapolation of cdf 
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« 


up 


S . 3 E xample 2: A Local Strain Analysis Problem 

in a Reliability Format 

A fatigue life estimate is required for a turbine disk. It is assumed that 
the fatigue crack will initiate at a notch in the rim provided for a blade as 
indicated in Figure 3.2. Stress and thermal analyses based on the definition of 
a duty cycle (one flight) will provide the normal stress S as well as the equi- 
valent stress v-g = K^S and the temperature ?t the notch as a function of time. 

A typical stress and temperature history shown in Figure 8.3 shows that the peaK 
stress will not, in general, occur at the same time as the peak temoerature. 

Generally speaking, in the event of cyclic plasticity at a notch, the local 
strain approach to fatigue life prediction (1, 2, 3) would be the method of pref- 
erence. The most serious limitation to such an approach in dealing with high 
temperature is that the strain life relationship must be based on a given constant 
temperature. 

Application of the local strain approach to high temperature fatigue prob- 
lems has considerable appeal. This approach to fatigue life prediction is now 
commonly used by mechanical designers in many fields. Analysis methods are well 
developed, including life prediction of components subjected to random loads 
and having notches which experience cyclic plasticity. While this method cannot 
explicitly treat the synergistic effects of creep and fatigue, life estimates 
may be reasonable for those components where creep effects may be small in com- 
parison to fatigue. This may be the case for example for some turbine hot section 
components where cycle lives of the order of 10“ a'^e anticipated. 

Following will be a demonstration of how reliability analysis methods, pre- 
viously described, can be applied. First, it is necessary to identify sources 
of uncertainty in the life prediction process. 

1. The use of the "duty cycle" to represent operational conditions. 

2. The temperature of the hot gases. 

3. The use of an equivalent temperature. 

Sehaviov 

1. Effect of strain hardening (or softening). 

2. Accuracy in the cyclic stress strain curve used. 


Notch in RiiTi or Turbine Disk Where Fatigue Crack 
is Assumed to Initiate 


Section of Turo’ne Disk 



Start Here 



Stress ill 




3. The use of Meuber's rule to obtain tne nominal stress- no t-.h st'^ain 

curve. 

4. Basic fatigue behavior as evidence by statistical scatter in strain- 
life data. 

5. Effect of surface finish, size effect, grain directions, etc. on the 
strain-1 ife curve. 

6. The assumption to include (or not include) mean stress relaxation. 

1. Accuracy of heat transfer analysis. 

2. Stress analysis procedures (including thermal stresses) to determine 
inelastic notch strains or elastic stress concentration factors. 

3. Procedure used to establish the stress-strain history at the notch. 

It will be assumed in this demonstration that the service nominal stress 

varies from zero to peak S as shown in Figure 8.4. The equivalent stress at the 
notch also varies from zero to a peak of = K^S where is the theoretical 
stress concentration factor. The equivalent and assumed constant temperature is 
denoted as T. 

Knowing the cyclic stress strain curve and applying Neuber's rule, the actual 
stress strain (a - c) curve at the notch can be established using standard methods 
of local strain analysis (1, 2). As illustrated in Figure 8.5, the resulting a-c 
history so established will define the strain range Ae and mean stress a^, the 
information needed to compare with the strain-life relationship. 


The strength of the material at a given cycle life N is defined by the general 
expression 




^ (2N) + e^(2N) 


(3.3) 
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in which is strain amplitude and S^, b, and c are empirical constants. The 
term in brackets is the familiar strain-life relationship (4.5). B is a random 
variable which accounts for the scatter in fatigue data. Using this form, it is 
assumed that the standard deviation of log (strain) given N is a constant (See 
Section 5.3 for further detail). 

The mean stress is a and will be a function of both temperature T and the 

m 

equivalent stress at the notch = K^S 

= a„(T,K.,S) (8.4) 

m m t 


Similarly the service strain amplitude will be a function of T, and S. 

= Ac/2 = C5(T,K^.S) (8.5) 


Consider operation to a given design life Failure is said to occur if 
"strength" is less than "stress". The event of failure is or 



[^f ■ 

L : 


(2No>'= 



< e5(T,K^,S) 


( 8 . 6 ) 


In this example, it is assumed that both and S are random variables re- 
flecting uncertainty in analysis procedures. Both and S will be the same for 
each cycle, but there is uncertainty in the value. Temperature T may be the most 
important variable but it is assumed to be constant just for the purpose of this 
demonstration. No theoretical difficulties would be introduced by considering T 
as a random variable. 

Data for this problem is summarized in Table 8.3. The stress strain history 
shown in Figure 8.6 obtained using Neuber's rule (see Reference 3) shows that 
only elastic strains are present after the first cycle. Assuming no mean stress 
relaxation, the stress strain history will cycle between points B and C as the 
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Table 8.3 


Random Variable 


Data for Example Problem 


Median Coefficient of Variation 


3 0.06 

40 ksi (276 MPa) 0.10 

120 ksi (827 MPa) 0.117 

1 0.18 


Constants 


fio 5000 cycles 

210 ksi (1448 MPa) 
b -0.07 

4 1-20 

c -0.70 

E 30,000 ksi (207 GPa) 


From basic properties of lognormal variates 


: = K.S 

e t 


cov(ag) =y 


(1 + C^Hl + C^) - 1 
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Figure o.5 

Determination of Actual Stress-Strain History at Notcn 
for Nominal Stress, S = 42 ksi 
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nominal stress cycles befv^len 0 and S. From a graphical analysis, ^See Figure 8.6 
for various values of S) it is shown that 


= ^-zl2 = K^S/2E = K^S/60,000 


(8.7) 


= 52.0 ksi (359 MPa) for any S 

Two methods of reliability analysis are used. In both it is assumed that K^, 
S and B have lognormal distributions. 

1. The Monte Carlo method; 

Fatigue strength (strain at failure) is defined by Equation 8.3. 

Noting that is constant, it follows upon substitution of the values of 
Table 8.3 that 

e, = 8(0.00467) 

Q 

Service strain is given by Equation 8.7. Failure occurs when 


8(0.00467) < K^S/60,000 

Isolating the random variables on the left hand side, the condition of 
failure is 

V = B/K^S < 3.6 X 10"^ 

u * 


Values of K^, S and B are selected at random from lognormal distributions 
having the parameters as given in Table 8.3. Then a sample of V is computed 
and its empirical distribution function is established as shown in Figure 
8.7. The figure shows the estimate of the safety index as 4.10 and the 

_5 

probability of failure as p^ = 2.09 x 10 . 
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INVERSE CDF OF F (.g) 


I 

i 


j- .9995 Figure b.7 
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2. The aenera 1 i red safety index; 

HASL was used to compute the generalized safety index, 3 = 3.97. A log 
transformation was made on all variables (See Section 7.4). The corresponding 
= 3.51 X lO"^. The output is shown in Table 8.4. This problem satisfies 
the assumptions of the lognormal format (See Section 7.2) and therefore HASL 
will provide the same result as the closed form expressions (Equations 7.20, 
7.21, 7.13 and 7.14). Moreover, the Monte Carlo solution should agree. The 
slight difference in results illustrates the error that is possible in extra- 
polating the distribution function into the tail regions. 

In both cases an index of structural performance was computed. The decision 
with regard to safety must now be made by the designer. 

As a footnote to this problem, consider the overall reliability of the disk. 

-5 

The estimate of fatigue crack initiation at a given site was p^ ~ 3.51 x 10 

But assume that there are n = 36 blades (i.e. n = 36 initiation sites). The key 

question is what is the probability of crack initiation at any site? 

This is a standard reliability problem, i.e. a series system where failure 

of the disk is assumed to occur if failure occurs in any component. System 

reliability is n n 

= li R. = n (1 - p.)" (8.8) 

s i=l ' i=l ' 

where R. = reliability of i^*^ component and p^. = probability of failure of i 
component. If p^ = constant, as it is in this example 

Rg = (1 - p)" (8.9) 

-5 36 

In this example R^ = (1 - 3.51 x 10 ) 

= 0.99874 
= 0.999 
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Table 8.4 

Output of HASL^^\ A Program for Computing the Generalized 

Safety Index for Example 2 


DESIGN VARIABLES 


VAPI ABLE 

(b) 

transformation' 

ME AN/HEDI4N 

(c) 

STD OEV/COV 

B 

1 


.lOOOOE«>Ol 

•ISOOOE^OO 

KT 

1 


. 30000E+01 

.60000E*01 

S 

1 


•60003EO02 

•lOOOOE^OO 



DESIGN POINT 


VARIABLE 

REDUCED VALUE 

BASIC VALUE 

B 

XP (1J» 

-3.33A5P 

X(l)> 

• 95132E«>00 

KT 

XR(2)» 

1. 1191^ 

X(2)« 

•320B2E^0l 

S 

XR (3)» 

1.86133 

X(3>« 

.A8152E+02 

SAPETY INDEX, 

BETA • 3* 

97860 



NOTES: 






(a) HASL is a digital computer program developed by P. H. Wirsching, University 
of Arizona, for computing the generalized safety index 

(b) A "1" indicates a log transformation on each variable 

(c) With a log transformation, the listed values are the Median and the Cov 
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Thus we would predict that only about 1 in 1000 of these disks would initiate 
a fatigue crack at the notch. 

3.4 Example 3: A Demonstration of the Efficiency of the Rackwitz Algorithm 

The Rackwitz iteration method, as described in Section 7.5, was oresented 
as a method for efficiently computing estimates of probabilities of failure 
for structural and mechanical components. The statistical distribution of 
each random variable is specified. The limit state can take any form, 
although risk estimates are expected to be slightly in error for non-linear 
cases. 

Following is an example which illustrates the efficiency of the Rackwitz 
algorithm as compared to Monte Carlo. In this example, computer costs for 
the Monte Carlo solution were six times greater. 

Consider a fatigue problem in which the fatigue strength is defined by 
the general strain-life relationship. Cyclic strain at a notch has a 
constant amplitude (not random), but there is uncertainty regarding its mag- 
nitude. Mean stress also remains constant, but its value is uncertain. 

m 

The limit state is written as 

I 

S " (8.10) 

The right hand side of the equation is fatigue strength. All terms have 
been defined previously in Section 5.3. Failure occurs if the strain (left 
hand side of Equation 8.10) exceeds the strength. 

I 

Values of the design factors are given in Table 8.5. Note that c , s^, 

I 

and are considered to be random variables. They have different 
distributions. In general it would be very difficult to compute the proba- 
bility of failure p^ using the exact form of Equation 7.5. However two 
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Table 8.5 


Values of Design Factors for Fxample Problem 

4 

Design Life, N = 10 cycles 
Modulus of Elasticity, E = 29,500 ksi 
Fatigue Strength Exponent, b = -0.108 
Fatigue Ductility Exponent, c = -0.540 


Random Design Factors 

Distribution 

Mean 

Standard Deviation 

1 

Fatigue Strength Coefficient, s^ 

* 

Lognormal 

148. 

.042 

1 

Fatigue Ductility Coefficient, 

* 

Lognormal 

.491 

.263 

Strain Amplitude, 

EVD 

.002 

.0003 

Mean Stress, a 

Weibull 

20. 

2. 


Values of the mean and standard deviation stated for the lognormal are really 
median and coefficient of variation respectively. 
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approximate methods will be employed to estimate p^. 

The Monte Carlo solution is presented on Figure 8.8. The safety index 
determined graphically is approximately 8 = 2.95. The corresponding p^ = 

0.0016. 

Solutions using the Rackwitz algorithm using SUPER-HASL is summarized 
in Table 3.5. Note that the results are essentially the same. In this 
example, the approximate costs of running the Monte Carlo program and SUPER- 
HASL were $6.00 and $0.70 respectively. Moreover, computer core requirements 
for SUPER-HASL are very modest. The program can be operated on a small 

computer. 

The quality and efficiency of the Rackwitz algorithm will be the subject 
of a continuing study. 


8.5 Summary Comments 

The three examples provided in this chapter demonstrate application of 
fatigue reliability methods. The approach chosen however depencs upon the 
statistical information available as well as the form of the limit state. 

This report is considered to be preliminary. Additional examples are 
being ronsidered and will be included in future reports. 
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SUPER-HASL* 


EFFICIENT ITERATIVE SCHEME FOR COMPUTING PROBABILITIES OF FAILURE^ 
R. RACKWITZ. Technische untversltiit, Munchen, Germany i 


INPUT 


DESIGN VARIABLES 




mean/meoian 

STD DEV/COV 

VARIABLE 




ES(E ) 

EVD 

.20000E-02 

.30000E-03 

s 

SF(SL) 

LOG 

.14800E+03 

.42000E-01 

SM(a„) 

WEIBULL 

. 20000E+02 

.20000E+01 

' m 

EF(ep 

LOG 

.49120E+00 

.26300E+00 


RESULTS 

DESIGN POINT 


VARIABLE 

reduced value 

BASIC VALUE 

ES 

XR(1) = 2.47612 

X(l) » .30365E-02 

SF 

XR(2) * -.26978 

X(2) = .14633E+03 

SM 

XR(3) * .08325 

X(3) = .20402E+02 

EF 

XR(4) * -1 .53516 

X(4) = .33024E+00 


SAFETY INDEX, BETA = 2.92711 
PROBABILITY OF FAILURE = .001716 


* j Kw fnlleen F Kelly, graduate student in Aero- 

Computer program developed ‘nn^yersity of Arizona 

space and Mechanical Engineering at the university 
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INVERSE COF OP F (~p) 



f 1 1 I I I 1 I 

> 1.25 0 1.25 2.50 3.75 5.00 6.25 


LIMIT STATE, X >*10'’ 
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APPENDIX 1 THE LOGNORMAL DISTRIBUTION AND PROPERTIES OF LOGNORMAL 
VARIABLES 

The lognormal format is proposed herein as a structural reliability 
model. In this regard Appendix 1 provides, for reference purposes, 
a detailed description of the lognormal distribution as well as certain 
properties of lognormal variables that are used for the design equatiors. 
Additional information on the lognormal is given by Benjamin and Cornell 
and Ang and Tang*. 

Given a random variable X. If Y » In X (or log^^ X) has a normal 
distribution, then X is said to have a "lognormal" distribution. The 
mean and standard deviation of X and Y are and (uy,ay) 

respectively. The probability density function (pdf) of X is 

f^(x) * exp (In X - uy)^ (Al.l) 

Cy X _ 2(?y 

for X > 0 
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The "statistical parameters" for the lognormal distribution are (uy^Oy) 
The mean and standard deviation of the lognormal variables are how- 

ever the median of X, denoted as X and the coefficient of variation of X, 
are often used. Relationships between parameters and moments are 
1 . Base e 


= In j In (1 + or Uy = In X (A1.2) 

Jy^ » In (1 + Cy^) (A1.3) 


* exp [uy + 2 <^Y J 

- 1 ) 


Cv » / e ' 


2. Base 10 


(A1.4) 

(A1.5) 


^°^10 '‘X ” 2 '"X 

ay- = 0.434 log^Q (1 + 

(uy ^ j log,Q (a//0.434); 
= 10 ^ 


/ (a//. 4347 
C, = /lO 


- 1 


1 


no 


A 


(A1.6) 

(A1.7) 

(A1.8) 

(A1.9) 


In general, it is more meaningful to use the median (rather than 
the mean) as a measure of central tendency for any random variable 
having a laroe variance. The median of X, denoted as X, is defined 
by the relation = 0.5C, where is the distribution function. 


A- 2 


Thus X defines the point below which 5C" of the population is expected 
to fall . 

Tifie expression To-' i.i Equations and » 1.6 can be derlted from the 
definition of the median 

0.50 = P(X ^ X) = P(ln X In X) 

= P(Y ^ In X) 

= Fy(ln X) (Al.lO) 

But, Y, the median of Y, is 

0.50 = Fy(Y) (A1.11) 

and it follows upon comparing Equations Al.lO and Al.ll that 


Recall that Y has a normal distribution and as a result of symmetry, 
the mean of Y equals the median, Uy = Y. Thus 


Uy = In X 


X = exp (uy) 


Equations A1.2 can be solved for 


exp iuy) >'1 ^ 


(A1.13) 


The mean of X can be expressed as a function of the median. Equation 
A1.13 can be written as, 


(A1.14) 


(A 1.1 5) 
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Dividing Equation A1.15 by A1.14, 


This relationship between X and proves to be useful because information 
on and is often available. But when probability computations 
are made, it is generally more convenient to use X. 

How to make probability calculations. 


P(X < x^) = P(ln X < In X ) 

y — 0 ' 


= P(Y ^ In x^) 


(AI.17) 


But Y is normally distributed with mean and standard deviation ( 1 ^^, Oy). 
Thus, 


Y -u. 

P(X < X ) = P(-^ 

0 Gw 


% • ‘•y 


= $ 



(AI.18) 


where $ is the standard normal distribution function. 

Replacing Uy with X using Equation 4.13, an expression which is generally 
more convenient to use follows. 


P(X 1 Xq) = $ 


In (x^/X) ] 



(A1.19) 


A similar expression for probabilities is valid when base 10 logs are 
used. 
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How to estijiate parameters from data. Given a random sample, of size 


n, of X, 


X = (X,,X-,..., XJ 


the data is converted to a random sample of Y. 


Y = (Yj.Y^..... V^) 


where = In X^. 

Then the estimates of Uy and Oy are* 




2 c 2 - 1 ? rv yi2 


(Ai.20) 


(A1.21) 


(A1.22) 
(A1 .23) 


Property of lognormal variates useful for design . Consider an 
expression involving several random variables, 

„ AX^Y^^ 

nd 


(Ai .24) 


where X, Y, L) are independent lognormal variates, b, c, and d are 
constants. It will be shown that V also has a lognormal distribution. 

A form for probability calculations on V will be derived. Taking the 
log of both sides, 

log V = log A + b log X + c log Y - d log U (A 1.25) 


Maximum likelihood estimators 
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Let, 


V = logv , X ' = log X, . . . etc 


(A1 .26) 


Then, 


V = log A + bX' + cY‘ - dU ' (A1 .27) 

Note that X', Y' and U' are normal. It follows from the addition 
property of normal variates that V is also normal. Therefore V 
has a lognormal distribution. 

To make probability calculations, it is necessary to find 
parameters (V, cTyi), Taking the expected value of both sides of 
Equation A1.27, 

Vy. * log A + buj^, + cuy, - duy, (A1 .28) 

But recall (Equation A1.13) that, = log V, ... etc. Substituting for 
each u. Equation A1.28 becomes 

log y = log A + log + log Y*^ + log (A1 .29) 


Taking both sides to the "e" power, 


i) = 


AX*^Y*^ 


(A1 .30) 


Taking the variance of both sides of Equation A1.27, 

2 . . 2 2 . ^2 2 ^ .2 2 /at -ji i 

ay, = b + C ly' -31 ) 

2 2 

But recall from Equation A1.3 that = In (1 + ) ... etc. It follows 

that 


A-6 


(AT .32) 


= In 0 + ^ - C,/)'^^ 


Then to compute probabilities, consider Equation A1.19, 


P(V £ v^) 


= 


I 

L -V' J 


(A 1.33) 


Example of the "Lognormal Format " 

Consider a design problem in which time failure of a com- 
ponent, denoted as T, is a random variable. T is known to be a 
function of random design factors S and (A description of the 

physical problem and what these terms mean is not important in 
this example.) 

T = 3^/.'P 

where a and b are constants. It is assumed that the median and 
COV of ; and A are known. The design life of the component is T^. 

The event of faHure of the component is (T < T^), and the 
probability of failure is. 


Pf ■ 1 T„) 

At this point, if we can assume that 3 and A nave lognormal distributions, 
then T has a lognormal distribution, and from £q A1.19, (note the reversal 
of T^ and t) 
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where T' = In T. Define the safety index 3 as. 


3 = In (T/T^)/o^. 


Then, 


= "(- 2 ) 

3 plays the same role in probabilistic design as the factor of safety does 
in conventional design. To compute 3, the parameters must be evaluated 
using a form of Equations A1.30 and A1.32, 

Of = “in T ° 





